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PART 1 


INTRODUCTIOIT 




lUS PEEEELEASE AEIffllMENT 


PREFACE 


On 15 MarcR 1911 , a splinter meeting of the Shuttle users ICD meeting was 
held during which attention was directed toward the errors involved in trans- 
ferring the Orhiter IMU alignment to the lUS guidance system. Questions arose 
regarding what errors were pertinent, their nature, and the alignment transfer 
accuracy achievable. Failure to align with sufficient accuracy apparently 
implied a need to install a star tracker on the lUS. 

Boeing had assumed a per-axis alignment transfer accuracy of 6.3 min (3 “’ ) 
in connection with navigation error analysis of four different lUS reference 
missions (reference l). This value was based on the understanding that the 
Orbiter IMU per axis alignment errors would not exceed 6.0 mn (3^) at the time 
of alignment transfer. After it was purported at the 15 March splinter meeting 
that the Orbiter alignment error might significantly exceed 6.0 mn, with the 
implication that the alignment transfer error would significantly exceed 6.3 mm, 
Boeing stated that a star tracker would be required on the lUS in order to 
achieve the reference mission requirements (reference 2). Thus, the 6.3 min 
alignment transfer accuracy appears to stand as the lUS IMU alignment accuracy 
requirement . 

NASA/JSC took the action to evaluate Orbiter/lUS alignment transfer. The 
first document (reference 3) under this action, titled "Orbiter In-Orbit 
Alignment Accuracy", dated 21 September 1911 , addressed the question of the 
Orbiter ‘s alignment accuracy, believed at the beginning of the task to be the 
major contributor to the overall alignment transfer error. The subject docu- 
ment, "lUS Prerelease Alignment", reports the results of analyzing alignment 
transfer accuracy. This second document shows that Orbiter in-orbit alignment 
accuracy is not a factor affecting lUS alignment accuracy, if certain procedures 
are followed. 

The basic analysis results are as follows. 

o Alignment of the Orbiter, per OFT procedures, followed by separate 

Orbiter /lUS alignment transfer procedures, just meets the lUS alignment 
accuracy requirement of 6.3 mu, if the elapsed time between Orbiter 
alignment start and alignment transfer completion is 20 minutes or less, 
the Orbiter alignment stars are essentially 90 degrees apart , and star 
images are restricted to the central 4 x degrees of the star tracker's 
field-of-view 

o The 6.3 min accuracy requirement is easily met by combining the Orbiter 

in-orbit alignment procedure, modified to remove sensor misalignment bias 
errors, and the Orbiter/IUS ali^mient transfer. In this case, expected 
lUS alignment accuracy is 1 6 min or better. 

The analysis results are more fully summarized in the next section. 
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SUMMARY OF RESULTS 


The analysis initially assumed that the Orbiter alignment and the Orbiter/IUS 
alignment transfer are performed as two separate procedures. It was then dis*- 
coTered that combining the procedures would lead to a substantial improvement in 
ms alignment accuracy. Results are summarized below first for the combined 
procedure and then for the separate procedures. 

Combined Orbiter Alignment and Orbiter/ITJS Alignment Transfer Procedure 

Combining the Orbiter in-orbit alignment procedure, modified to remove body- 
fixed sensor misalignments, and the Orbiter/IUS alignment transfer procedure, as 
prescribed below, leads to an lUS per-axis alignment accuracy of 1.6 min (3a) or 
better. The combined procedure is as follows: 

(1) Orbiter takes alignment sighting on star #1, using either of its two 
star trackers, recording star tracker and Orbiter IMU gimbal angle measurements. 

IDS attitude from the lUS strap-down IMd system is simultaneously recorded in 
IDS flight computer. 

(2) Orbiter rotates 180 degrees about star #1 line of sight (LOS) and then 
takes another alignment sighting (same star tracker) on star #1, again recording 
star tracker and Orbiter IMD gimbal angle measurements. lUS attitude is again 
simultaneously recorded in the IDS flight computer. 

(3) The first and second set of Orbiter measurements are averaged, removing 
the body-fixed sensor misalignment effects from star #1 measurements. In addition, 
the eigenvector (eigenvector #1) associated with the 180 degree rotation is com- 
puted in both the Orbiter and IDS flight conputers. The eigenvector essentially 
represents the axis of rotation. 

(4) The Orbiter selects alignment star #2 and repeats (1) , (2) , and (3) , 
using either of its two star trackers. This yields an averaged star measurement 
on star #2 (Orbiter computer) and eigenvector #2 (in Orbiter and IDS computers). 

(5) The Orbiter measurements, expressed in the Orbiter’ s IMU stable member 
inertial coordinate system, and the IDS measurements, expressed in the IDS inertial 
coordinate system (orientation unknown at this point), are jointly processed. 

(The implementation approach for computer processing of the measurements has not 
been definitely established at the present time, but it is understood that the 
Orbiter measurements will be supplied to the IDS flight computer where the align- 
ment computation will take place.) The end result of the processing is a 3 x 3 
Diatrix transformation that relates the IDS unknown inertial coordinate frame to 
the desired IDS inertial navigation coordinate frame (such as the M50 coordinate 
frame) . Applying the matrix transformation to IDS body attitude (one shot compu- 
tation) constitutes the IDS alignment. 

If the alignment stars are 90 degrees apart, the per-axis error of the IDS 
alignment is 1.6 min (3a). If the Orbiter star tracker measurements are restricted 
to the central 4x4 degree field of view (full field of view is 10 x 10 degrees) , 
then the per-axis error is 1.0 mh (3a). 
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Table 1. lUS Alignment Accuracy, Combined Procedure 


Star Tracker FOV 

lUS Alignment Accuracy ( 3 cr) 

10 X 10 deg 

1,6 min 

4 X 4 deg 

1.0 min 


Note; Alignment stars are 90 degrees apart 

The lUS alignment error is due to (l) the Orbiter star tracker and (2) the 
ITJS IMU gyros . Orbiter IMU readout , drift , and alignment errors essentially do 
not impact the lUS alignment accuracy, given the procedures outlined above. 

The average per-axis alignment error degrades by the factor K = (l + 2csc A)"^/ 
■^3 j Where A is the subtended angle between the alignment stars. When A = 90 degree' 
K = 1. For 60 degrees <A<120 degrees, IKl.l. 

Separate Orbiter Alignment and Orhiter/IUS Alignment Transfer Procedures 

When the Orbiter alignment and the Orbiter/IUS alignment transfer are accom- 
plished via separate maneuvers and procedures, then the Orbiter IMU readout, drift, 
and alignment errors directly impact the lUS alignment acciiracy. Two cases were 
analyzed with the following results. 

Case A 


The Orbiter is aligned per OFT in-orbit alignment procedure (reference 4), 
thus the error effects of Orbiter body-fixed sensor misalignments are not removed. 
Alignment is subsequently transferred to the lUS via rotations about two axes 90 
degrees apart in inertial space. lUS alignment accuracy is presented below for 
(l) rotation magnitudes of 90 and l80 degrees and (2) Orbiter star measurements 
over the full 10 x 10 degree tracker field of view (FOV) and restricted to the 
central 4x4 degree FOV. The effects of Orbiter gyro drift (.035 deg/hr, lo") , 
which depend on elapsed time after Orbiter alignment, are also indicated. 
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TalDle 2. Case' A lUS Alignment, Orbit er Sensor Body-Fixed Misalignment 
Errors Not RemoTed* 


Elapsed Time* 

Alignment Transfer 

lUS Per-Axis Alignment 


Rotation Magnitude 

Accuracy (3 o') 



10 X 10 deg FOV 

4x4 FOV 

0 min 

90 deg 
l80 deg 

6.3 mih 
6.1 min 

6.0 ^iS 
5.8 min 

20 min 

90 deg 

6.7 min 

6.4 min 


l80 deg 

6.5 min 

6.2 min 

itO min 

90 deg 

7-6 min 

7.3 mu 


l80 deg 


7.2 min 

60 min 

90 deg 

8 . 9 min 

8.7 min 


180 deg 

8 8 min 

8.6 mn 


Note 

o Elapsed time is period between start of Orbiter alignment and end of 
Orbiter/lUS alignment transfer 

o Angle between Orbiter alignment stars assxaaed 90 degrees. 


OJilGlNAl^.l^Gte IS 
OF POOIt QUAinPY 

' ti 
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Case B 


The orbiter is aligned per OFT procedure modified to remove sensor body-fixed 
misalignment bias errors from the Orbiter star measurements. (This is the pro- 
cedure described in step ( 2 ) of the combined alignment transfer procedure described 
earlier.) Alignment is subsequently transferred to the ITJS via rotation about 
two sixes 90 degrees apart in inertial space. The lUS alignment accuracy is pre- 
sented below . 


Table 3 - Case B lUS Alignment, Orbiter Sensor Body-Fixed 
Misalignment Errors Removed* 


Elapsed Time* 

Alignment Transfer 
Rotation Magnitude 

lUS Per-Axis Alignment 
Accuracy (3 o') 

10 X 10 deg FOV 

4 X 4 deg FOV 

0 min 

90 deg 

3 -T min 

3.4 min 


180 deg 

3.1 mn 

2.8 min 

20 min 

90 deg 

4.3 min 

4.0 mn 


180 deg 

3 -T min 

3.5 

^iO min 

90 deg 

5.6 mn 

5.4 ^ 


180 deg 

5.2 min 

5 0 min 



^ 


60 min 

90 deg 

7.3 min 

7.2 min 


180 deg 

7.0 min 

6.9 min 


*Note * 

o Elapsed time is as defined for Table 2 

o Angle between Orbiter alignment stars assumed 90 degrees 

i ' ' t ' . 

For Case B'the alignment transfer rotations are maneuvers separate from the 
180 degree bias removal rotations about the Orbiter alignment star LOS's. It 
would seem that if the bias removal procedure is adopted, there would be little 
reason not to combine the Orbiter alignment and Orbiter/lUS alignment transfer 
procedures . Cobibining the procedures reduces the amount of maneuvering time re- 
quired, and it improves the lUS alignment accuracy to that reported in Table 1 . 
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ALIGNMENT TEANSFER: THE BASIC IDEA 


Fox many lUS missions, the Orhiter will transport the lUS into orbit in a 
powered down state. One of the necessary steps in preparing the lUS for release 
is aligning the lUS’s strapdown inertial guidance unit (IMH) to the basic inertial 
coordinate frame of reference chosen for the mission. Targeting data and the 
gravity model are stored in the lUS flight computer in such a basic reference 
frame. We can assume this inertial frame to be the same as the Orbiter's, i.e., 
the M50 inertial frame, without loss of generality. 

At some point prior to lUS release from the Orbiter, the lUS’s IlflJ and flight 
computers are powered up. At this time, the flight computer has no idea of the 
lUS’s orientation relative to the M50 coordinate frame. The lUS flight computer 
does begin, however, to track lUS attitude relative to the inertial attitude 
existing at the Instant the attitude computations were initiated. Thus, the lUS 
has an inertial reference frame, but it doesn't know the frame orientation rela- 
tive to the M50 frame. The lUS inertial frame, at this point, is unknown. 

Aligning a strapdown IMU consists of determining the orientation of the lUS’s 
unknown inertial frame relative to the M50 frame. Then, the accelerations sensed 
along the body axes by the IIJS's IMG can be expressed in M50 coordinates and 
combined with the gravity model accelerations as the basic inputs to the navigation 
computations . The alignment itself is expressed mathematically in the lUS flight 
computer as a 3 x 3 matrix transformation. The alignment procedure addresses the 
problem of determining the 3x3 matrix transformation via appropriate sensor 
measurements and vehicle maneuvers . 

The basic idea behind transferring Orbiter alignment to the lUS is the 
following. The Orbiter, with the lUS firmly attached in the payload bay, performs 
rotations about two different spatial axes. The rotations are jointly sensed by 
the Orbiter and the lUS, affording two common lines of reference in inertial 
space. For the Orbiter, the two reference directions are expressed in the Orbiter's 
inertial coordinate frame. For the lUS, the same two reference directions are 
expressed in the lUS's ^3nknown inertial frame. Since the reference directions are 
common to both the Orbiter and the lUS , it becomes a simple matter to compute 
the orientation of the ITJS^ unknown frame with respect to Orbiter's frame. 

If the Orbiter alignment and the Orbiter/IUS transfer maneuvers take place 
separately, then the alignment transfer error will be the sum of the Orbiter 
alignment, Orbiter IMU, and the lUS IMU errors. It will be shown that the Orbiter 
in-orbit alignment procedure and the Orbiter/lUS transfer alignment procedure can 
be combined with the consequence that only the Orbiter star tracker errors and 
the lUS IMU errors impact the lUS ali^ment. For this latter situation, it turns 
out that lUS alignment is substantially more accurate than the Orbiter's alignment. 

Analysis Approach 

At the present time, Orbiter/lUS alignment transfer procedures and calculations 
have not been explicitly baselined. Various approaches appear feasible. For this 
analysis, the liberty was taken to adopt a simple deterministic computational 
approach. The equations involved are derived in the text of this report. Such an 
approach is easy to understand, provides a framework upon which the error analysis 
can be performed, and furnished an accuracy benchmark against which other approaches 
can be conpared. 
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The hasic error sources considered in the analysis are the sensor error 
soui^ces, discussed in a later section. These are the errors associated with the 
Orbiter IMD's and star trackers and the lUS’s strapdoxra IMU, Two other potential 
error sources exist, but at the present time quantitative data is not avaxlable 
to evaluate their significance. Thus, they were not included in the analysis. 
These error sources are: 

o Data processing system implementation (mechanization) errors . The 

principal error source here would be timing errors associated with time 
tagging of measurement data. This error source is not significant, if 
the Orbiter and lUS master time references are known relative to each 
other within several milliseconds. 

o Variations of IDS orientation relative to the Orbiter, when the lUS is 
attached to the Orbiter payload bay. Accurate Orbiter/IUS alignment 
transfer is predicated upon the assunqjtion that the Orbiter and IDS 
rotate as a single unit during alignment transfer maneuvers. Strictly 
speaking, changes in the lUS’s navigation base relative to the Orbiter’ s 
navigation base, from measurement to measurement, will introduce errors 
into the alignment transfer process. For example, a .1 degree relative 
shift in IDS navigation base orientation, say due to body flexing or 
Orbiter/IDS attach points that are not rigid, might yield a 6 inin error 
in the IDS alignment. 

It would seem very important to verify, or be quite confident, that the 
IDS relative orientation does not change significantly from measurement to 
measurement . 

This document will first address IDS alignment transfer via separate 
Orbiter alignment and Orbiter/IUS alignment transfer procedures. Building 
on the analysis results and procedures developed thereby, the analysis will 
then address IDS direct alignment via combined procedures . 

Comment on the Appendices 

Appendix A is a detailed discussion of eigenvector calculations pertaining 
to IDS alignment. The other appendices (B,C,D, and E) are sections from the 
first report under this analysis task, "in-Orbit Alignment Accuracy", reference 
3. They are Included for completeness and convenience of the reader, since the 
subject report makes a number of references to the first report. 
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NAME* Aries mean of 1950, CartC'Sian, coordinate system 

ORIGIN The center of the earth 

ORIENTATION: The epoch is the beginning of Besselian Year 1950 or 

Julian ephemeris date 2433282.423357 

The Xj^j-Yj^ plane is the mean earth's equator of epoch. 

' The Xf^ axis is directed towards the mean vernal equinox 

of epoch. 

The axis is directed along the earth's mean rotational 
axis of epoch and is positive north. 

The Yj^ axis completes a right-handed system. 

CHARACTERISTICS. Inertial, right-handed Cartesian system. 


Figure 4. 2.1-7. Aries Mean of 1950 Cartesian Coordinate System 


4-2-9 ' 

sn7r>-5U-ooi3 

17 0*70 
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NOTES ■ 


Ya 


Xg = X GYRO 


= 

^G = 


I. A. 

I. A, J 


VERTICAL GYRO 




Y GYRO 
Z GYRO I. A. 

RED GYRO I. A J 
Z ACCEL SENS AXIS 
X ACCEL SENS. AXIS 

Y ACCEL SENS. AXIS 


AZIMUTH GYRO 


) 


*'G 

= 

= 

- GIMBAL TORQUER 
= GIMBAL ANGLE RESOLVER 


SINGLE AXIS ACCEL 


DUAL AXIS ACCEL. 




Figure 4. 2. 1-4. Stable Member (IMU Cluster) Coordinate System (Sheet 1 of 2) 
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NAME 


STABLE MEMBER (IMU) 


ORIGIN THE INTERSECTION OF THE INNERMOST GIHBAl AXIS AND THE 

MEASUREMENT PLANE OF THE XY TVO AXIS ACCELEROMETER 

ORIENTATION THE AXIS IS COINCIDENT WITH THE INNERMOST GIHBAL AXIS 

THE Xc AXIS IS DETERMINED BY THE PROJECTION OF THE X 
ACCELEROMETER INPUT AXIS (I A) ONTO A PLANE ORTHOGONAL TO 
Zj. Yj. COMPLETES A RIGHT-HANDED TRIAD. 

IN A PERFECT IMU. WITH ALL MISALIGNMENTS ZERO, THESE 
RELATIONSHIPS HOLD 

THE X ACCELEROMETER AND X GYRO IAS ARE PARALLEL TO THE Xj,AXlS 

THE Y ACCELEROMETER AND Y GYRO IAS ARE PARALLEL TO THE Y-AXIS 

the 2 ACCELEROMETER AND Z GYRO IAS ARE PARALLEL TO THE 2^ AX IS 

CHARACTERISTICS NONROTATING, RIGHT-HANDED, CARTESIAN SYSTEM 

THE REFERENCE ALIGNMENT FOR THE GIMBAL CASE SHALL BE 
defined WITH THE FOUR GIMBAL ANGLES AT ZERO AND WITH 
THE VEHICLE IN A HORIZONTAL POSITION. IN A PERFECT 
IMU, WITH ALL MISALIGNMENTS ZERO AND WITH ALL GIMBAL 
ANGLES AT ZERO, THE FOLLOWING RELATIONSHIPS HOLD 

THE OUTER ROLL AXIS AND THE X^. AXIS WILL BE PARALLEL 

TO Xjjg. POSITIVE Xj. WILL BE IN THE FORWARD DIRECTION 

POSITIVE ROLL GIHBAL ANGLES WILL BE IN THE SENSE OF A 
RIGHT-HANDED ROTATION OF THE GIHBAL CASE RELATIVE TO 
THE PLATFORM ABOUT THE PLUS OUTER ROLL AXIS 

THE PITCH AXIS AND Y- WILL BE PARALLEL TO Y.,_ POSI- 

t. ND , 

TIVE Yj, WILL BE TO THE RIGHT OF AN OBSERVER LOOKING 

FORWARD IN THE VEHICLE. POSITIVE PITCH GIMBAL ANGLES 
WILL BE IN THE SENSE OF A RIGHT-HANDED ROTATION OF THE 
GIMBAL CASE RELATIVE TO THE PLATFORM ABOUT THE PLUS 
PITCH AXIS 

THE INNER ROLL AXIS WILL BE PARALLEL^ TO THE OUTER ROLL 
AXIS, WITH THE SENSE OF ROTATION THE SAME AS FOR THE 
OUTER ROLL AXIS 

THE AZIMUTH AXIS AND WILL BE PARALLEL TO 

POSITIVE WILL BE DOWN RELATIVE TO AN OBSERVER STANDING 

IN THE VEHICLE. POSITIVE AZIMUTH GIMBAL ANGLES WILL BE 
IN THE SENSE OF A RIGHT-HANOEO ROTATION OF THE GIMBAL 
CASE RELATIVE TO THE PLATFORM ABOUT THE PLUS AZIMUTH 
AXIS. 

*NB’ %B’ CARTESIAN COMPONENTS OF THE NAVIGATION 

BASE COORDINATE SYSTEM , 

Figure 4. 2,1-4. Stable Member (IMU Cluster) 

Coordinate System (Sheet 2 of 2) 
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Zb 


NAME: Navigation base system 

ORIGIN: At the mutual Intersection of: 

(A) A plane parallel to the orblter plane of symmetry, 

14 inches left of center 

*(B) Plane of top surfaces of mounting pads for IMU 1 (left) 

*(C) Plane of vertical surfaces of aft pads for IMU 1 

ORIENTATION: Yjjb lies along the intersection of planes (B) and (C) , 

positive out the orbiter right wing. 

XuB lies in plane (B) , perpendicular to Y^b and positive 
forward. 

ZjjB completes the right-handed system. * 

CHARACTERISTICS: Rotating, right-handed' Cartesian 

*As determined by an alignment fixture ORlGrlNAL PAGE 

<^F p<?oR QUALrrx 

Figure 4. 2. 1-3. Navigation Base Coordinate System ' ’ 
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COMPUTE RESOLVER READINGS IN RADIANS 
OR = (ORANG) (1/RADDEG) * 

P « (PANG) (1/RADDEG) 

AZ = (AZANG) (1/)IADDEG) 

IR = (I RANG) (1/raDDEG) 


* WHEN CALLED BY 
STAR TRACKER SOP 
ORAVG21, PAVG21, AZAVG21, AND 
IRAVG21 REPLACE ORANG, 

PANG, AZANG, AND I RANG 


FORM NAV BASE TO CLUSTER TRANSFORMATION. 

1 " o' 'l 0 o' CP 0 

^ SAZ CAZ 0 0 CIR-SIR 0 1 

_0 0 ij Lo SIR CIrJ -SP 0 




I - 


Figure 4 6. 2. 8-1. TNB_CL Flow Diagram 
4.6-127 
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4,66 STAR_rR^CKE5_SUBSYSrE^_0PERATING^gR0GRAM 

This section defines the detailed functional requirements and 
formulations for the Star Tracker (ST) Subsystem Operating Program 
(SOP) , The ST SOP defines software associated with ST moding, seif- 
test, failure identification, star tracking, target tracking and IMU- 
to-sr alignment, 

4’. 66.1 ST REaUIRE'lENTS OVERVIEW 

4.66.1.1 

The Orbiter ST is a strapped-dovm, wide field of view (FOVJ image- 
dissector electro optical tracking device. It is used to obtain 
precise angular measurements of selected stars and sun illuminated 
orbiting elements (targets) . 

Two ST's are physically located on the Orbiter navigation base. The 
-? axis ST centerline is approximately 10,5° from the Orbiter -? axis 
and the •z axis ST centerline is approximately 3° from the Orbiter -Z 
axis. Figure 4.66.1.1-1 depicts ST and IMU placement on the 
navigation base. The ST mounts on the underside of the navigation 
base while a light shade and viewing window are mounted on top of the 
navigation base. The ST interfaces wirh the GPC via the serial 
digital mulfciplex/demultiplex input/output data channel. 

4.66.1.1.1 ST^£e 2 .f ormance^C ha £act grist ics - 

The Orbiter ST»s exhibit the following performance characteristics and 
operating features : 

* "rhe ST'S will acquire and track the 153 brightest 
stars based upon the s-20 star intensity scale. The ST 
sensitivity threshold can be adjusted via GPC control. 

Ei§ld_of_View - The ST's total field of view is a minimum 10 
degree by 10 degree square. The ST's can also be commanded to 
search a 1 degree by 1 degree square field within the field of 
view. 

fi.g cu ;acy - The ST* s total error in measurement of star or target 
angles does not exceed 1 arc minute (1 sigma). Star or target 
magnitude measurement errors do not exceed an absolute maximum of 
- p ij6 , ^magnitude . 


4.66-1 
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Figure A. 66. 1.1-1. ST I>^U TIAV Base Orientation 
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ERROR MODELS 


The sensor errors pertinent to alignment transfer are presented in this 
section. Table 4 lists the Orblter BIU errors. These errors were discussed in 
some detail in the first report of this task (reference 3) . This discussion is 
repeated for the convenience of the reader in Appendix B to this report. 

The Orblter alignment error analyzed in reference 3 (also discussed in 
Appendices C and E) is presented in Table 5. Two situations are addressed; 

o Orblter is aligned in accord with OFT in-orbit alignment procedures, 
reference 4. 

o Orblter is aligned per OFT procedure modified to remove the effects of 
body fixed bias effects in the star tracker and UKJ measurements. This 
is accomplished for each of two star sightings by (1) taking star measure- 
ment, (2) rotating Orbiter 180 degrees around LOS to star, (3) taking the 
second star measurement, (4) averaging the two measurements to remove the 
bias effects. This modified alignment procedure is described and analyzed 
in reference 3 (also discussed in Appendix D) . 

It is seen that removing the bias effects from the star sightings materially 
reduces the alignment error. 

Table 6 presents the lUS IMU error sources significant to the Orbiter/lUS 
alignment transfer. This error model is based on information received from the 
Boeing Company. Note that Table 6 does not include such errors as "nav base to 
IM(J alignment error." This is not because the lUS IMU is perfectly aligned on 
its navigation base plate. It is because, rather, we are not depending a priori 
on any particular orientation of the IMU relative to its nav base when performing 
the Orbiter/IUS alignment transfer procedure. If the lUS IMU alignment were 
carried out using an lUS star tracker, then the precision of the mounting align- 
ments of both the lUS IMU and lUS star tracker would become significant. 

Redundancy 

The lUS IMU design contains redundant gyros (and accelerometers) . Such 
redundancy, when taken into account, should lead to a net reduction in IMU error. 
This analysis will not address redtindancy effects, since the lUS IMU turns out to 
be a minor error contributor to Orbiter/IUS alignment transfer. 
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TABLE 


ALIGW4ENT ERROR MODEL 

ERROR SOURCE SYMBOL ^ VALUE/AXIS (ig) 


OrLiter IMU 


Gyro Drift 

Nav Base Ref to Mounting 
Pads 

IMU Case to Pads 

Case to Outer Roll Gimbal 'SQr.^-o 

, OUK 


.035 deg/hr 

6o 

20 sec 

28 


Body 

Fixed 

Biases 


Non-Orthogonality 

Resolver °R 


Star Tracker 

Horizontal, Vertical 
Measurements 


30 sec 53 sec 
RG3 


h2 sec 


68 sec 
RSS 


Tracker to Nav Base Ref 


60 


NST 


60 sec 


Body Fixed Bias 


RSS (not including gyro draft) e 


11 (l a) 

5.T t^n (3 a) 


* Note These symhols are employed in the error analysis, reference 3. 
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Table 5. Orblter IMJ Alignment Error 


Alignment Error, per OFT Procedures — — <—114 sec (la) 

5.7 TOui (3a) 

. Alignment Error, per Modxfxed Procedures — '48 sec (la) 

2,4 MA (3a) 


Table 6. lUS IMJ Gyro Error Model (la), per Axis 


Uncompensated Random Drxft .009°/hr 

90 day bxas imcertaxnty ,007°/lir 
Thermal cycle stabilxty .005 /hr 
Shutdown repeatxbxlxty .002^/hr 
Contxnuous operating random .001 /hr 


Scale Factor Error- 


45 ppm 


90 day uncertaxnty 25 ppm 

Linearity 37 ppm 

Misalignment Stability 10 s^ 


Error Analysis Technique 

All the various error sources result in small angular errors, which are 
small rotations. As discussed in reference 3, small rotations can be expressed, 
to first order, as vectors. The vector magnitude is the angular error magnitude 
in radians. The vector orientation is the axis about which the angular error 
takes place. In general, individual angular error vectors have different mag- 
nitudes and orientations . The total error (vector) is the vector sum of the 
Individual error vectors . 

* r 

Each angular error vector can obvicvusly be 'resolved in X, Y, and Z components 
in any particular coordinate system of interdst. ^ ^ 

For a given error source, we will assume the per-axis (X, Y, and Z) error 
statistics to be the same and also uncorrelated from axis to axis. This seems 
to be a reasonable assumption, based on the data at hand. In general, the error 
ellipsoid corresponding to such an, error distribution is a sphere. Thus, the 
per-axis error statistics are invariant to coordinate' frame transformations., We 
can, therefore, assume that the per-axis error statistics presented in Tables 4, 

5, and 6 pertain to the same coordinate frame, with said frame being whatever 
coordinate frame we ^choose to perform the error analysis. 
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The star tracker is an exception to the above since xts "error ellipsoid" 
is really an ellipse in the plane perpendicular to the tracker boreaxght axxs. 
This will not present a problem in the analysis of per-'axis errors, since we will 
choose the analysis coordinate frame Z axis to be coincident with the tracker 
boresight axis. The tracker measurement errors (X and Y axis errors) are then 
directly additive to the other error source X and Y axis errors. Z axxs errors 
from any source (angular errors around the tracker boresight axis) to first 
order have little effect on alignment accuracy. 


ORIGINAL 
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PART 2 


ORBITER/IUS ALIGNMENT TRANSFER 


SEPARATE OKBITER ALIGNMENT 
AND ORBITER/IUS ALIGNMENT 
PROCEDURES 


original pagIjis 

OF POOR QUALITY 
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ElGEimCTOR INERTIAL REFEBMCE DIRECTION 


When the Orhxter/IUS performs an alignment transfer rotation, the angular 
velocity vector at any xnstant of txme represents a reference direction in 
inertxal space. If the angular velocity vector were sensed jointly by the Or- 
bxter and the lUS at the same time point, a common inertial reference direction 
would then be known. In practice, because of sensor errors and the nature of 
the sensor data (e.g., the Orbiter IMU provides angles, not rates), it would 
probably be necessary to collect IMU data over some period of time. The inertial 
reference direction would then be estimated from the data by some technique such 
as Kalman filtering. 

The analysis of this document does not attempt to work with angular rates. 

The basic measurement data is assumed to be body attitudes existxng (1) at the 
beginnxng of the rotatxon maneuver and (2) at the end of the rotation manuever. 
This attitude data is directly avaxlable to both the Orbxter and the lUS flight 
computers . 

There are an xnfinxte number of ways that the Orbxter/lUS could reorient from 
the inltxal attitude to the final attitude. In general, the instantaneous axis 
of rotation would vary throughout the maneuvering. However, there exxsts one 
"ideal" rotation of 180 degrees or less about fixed axis which would accomplish 
the given reorientation. This fixed axis, which can be easily calculated as a 
functxon of the Initial and final body attitudes, wxll be taken as the xnertial 
reference dxrection. This axis is independent of the actual maneuvering enqjloyed 
to reorient from the xnxtial attxtude to ending attitude. 

Mathematically, body attitude orientatxons are represented by orthogonal 
matrxx transformations . An orthogonal transformation has one independent eigen- 
vector whose directxon, xt turns out, corresponds to the fxxed "ideal" axis of 
rotatxon. Thus, the inertial reference direction xs represented by the eigen- 
vector of the 3x3 matrix transformation relating the initial and ending body 
attitudes . 

Rep re sent atx on of Body Attitude 

In this analysis vehicle body attitudes will be represented by 3 x 3 matrices 
of body axis dxrection cosines. Let A be such a matrix. 




^11 

t 

/ 

^12 1 

^13 

A 


^21 

1 

! 

®22 J 

^23 


r ^ 

_^31 

1 

H 

i 

^32 i 
1 

®33 


The three columns of A are three orthonormal unit vectors representxng respectively the 
roll, pitch, and yaw body axeSj_^left to right) as resolved in th|^gxven coordinate 
frame. Component a^j^^ is the x component (X, Yy or Z) of the unxt vector 

(roll, pitch or yaw) . It turns out that A xs an orthogonal matrix, with positive 
'determinant. 
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Let represent the mitral body orientatxon, expressed in a selected inertxal 

coordxnate frame, and represent the endmg body orxentatxon, expressed xn the 

same frame. We define a matrix transformation C relatxng A^ and A^. By defxnition 
we have; 




.-1 


C always exxsts, since A always exxsts; for an orthogonal matrxx, the inverse 
matrxx xs the transposed matrxx (A ^ = A^) . Thus ; 


C = A^A- 

C is seen to be an orthogonal transformatxon sxnce xt is the product of two orthogo- 
nal transformatxons . By defxnxtion of an exgenvector we have: 

Gd = d 


where d^ is the eigenvector of C. Likewxse it can be seen that: 


C^d = d 


Subtractxng yields: 

(C^ - C)d^ = 0 


Writing thxs equation xn component form yxelds: 


Hence : 
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We see that the equations above are satisfxed by; 


^ (632 



^2 ^ 


- k (Cji - 


where k is a constant. Since ^ has unit length, 

,2 . , ,2 


t *^32 ^^ 23 ^ ^^13 ^^ 31 ^ ^ ^*^21 “ ^ 12 ^ ] ^ 


Thus the eigenvector of C can be determined by picking the appropriate components 
out of the matrix C, coii 5 >uting k, and then forming 

We will see xn the error analysis that a mathematical singularity exists 
in the above solution for rotations of exactly 180 degrees . This does not preclude 
using this solution in the error analysis, which we wish to do because of its 
simplicity of form. In Appendix A, alternate methods for coD 5 )uting ^ are 
presented . 
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ALIGNMENT TRANSFER EQUATIONS. UNCOMBINED PROCEDURES 


We now develop the eqtiations fox determining lUS alignment wherein the Orbiter 
alignment and lUS alxgnment procedures are not comha^ned. 


Body Cosine Matrxces 


We express the Orbiter^ s attitude matrix A in Orbiter IMU stable member 
(present cluster) inertial coordxnates . 


A 




where 


T^ = orthogonal transformation (3x3), nav base reference 
to IMU outer roll gimbal. Known nominally, constant 
during rotation maneuvers. 

PC 

~ orthogonal transformation (3x3), outer roll gimbal 
to present cluster (stable member) frame. 

A = Orbiter body attitude matrix xn stable member frame. 


The equation above expresses the Orbiter ’s attitude in terms of IMU gxmbal 
PC 

angles, since T^, is a function of the gxmbal angle readings. 

KJj 

Corresponding to A is the lOS attitude matrix U, as expressed in the lonknown 
lUS inertial coordinate frame. 


First Rotation Maneuver 


The Orbiter/mS performs its first of two rotation maneuvers. In doing so, 
the Orbiter transitions from attitude A^^ to attitude A 2 . Likewise the lUS transi- 
tions from to U 2 . 

For the Orbiter, we compute the ideal rotation matrix C (discussed earlier) 
from A^ and A 2 with 


Since 


NB 


_ /^1 B1 
1 EL NB 

initial orientation 

PC. __ 

^2 EL NB 

final orientation 

C-AjAJ 

we have 

PC„ _ 

C = T ^ 

EL PC^ 

(stable member frame) 

Ls eliminated in 

the computations . Over 


ORIGINAL ^^aIjE IS 
OF POOR quality 


maneuvers . 


EL 

T.._ can be considered constant (reference 3) . This assiames the same 
NB 
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IliU provides tke angle measurements determinijig and 3 ..e., Orbxter redundancy 
management does not switch IMU’s -during a rotation maneuver, 


From Cj the eigenvector (corresponding to the first rotation axis) is 
computed using the appropriate equations presented later in this document. 

Using the Orhiter^s alignment matrix (determined by Orbiter alignment 

procedures), which transforms from the stable member frame to the M50 frame, 
is transformed to M50 coordinates . 



For the lUS, compute ideal rotation marrix F from .U^ end U 2 • 


From F, the eigenvector is computed, using the appropriate equations 
presented later in this document. and are nominally the same vector 

quantity in inertial space. However, they are expressed xn different coordinate 
frames . 

Second Rotation Maneuver 


The Orbiter/IUS performs the second rotation maneuver. Employing the same 
type calculations as above yields and the rotation eigenvector expressed 
in the M50 frame and the lUS unknown frame respectively. We assume and are 
not colinear (accordingly likewise and _g^) . 

Computation of lUS Alignment 

Form the following orthonormal triads; 


e^ = unit x ^ 2 j 

= unit X ^ 2 ] 

“2 “ * -1 

^2 ~ -^3 ^ -S-T 


Form the following matrices E and G from column vectors 
-S.2’ ft.3* 


® 1-^1 ^ 3 ] 

i .t'. 


^ % %] 


j The desired lUS alignment transformation matrix is which transforms any 

vector in the lUS unknown frame to the M50 frame. Clearly, transforms 

to to and to Hence: 


E = 


Solving for 3^*^ yields: 

=EG*= 
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lUS alignment taJces place t»y transforming thn IU3 body attitudes TJ from the 
lUg unknown frame to M50 coordinates. 

K = 


where K is the lUS body attitude matrix expressed in M50 coordinates. This is 

a one-time computation at some point in time subsequent to determining 

Once the transformation is made, the lUS flight con^uter will proceed to update 
K in the M50 frame, using its strapdown IMO inputs. 

As will be seen in the error analysis, the errors in per this iincombined 

procedure are due to the Orbiter alignment error in the errors in determ- 

ining and d^ due to IMO readout errors, and the errors in determining and 
£2 due to lUS IMO errors. 


Comment on Orbiter IMP Redundancy Management 

The Orbiter has three IMO’s. The on-board redundancy management function 
middle— value selects one IMO among the three IMO’s for input to the flight com- 
puters. As stated above, it would be important that the same IMO provide the 
measurements prior to and at the completion of a rotation maneuver. It is under- 
stood that the crew can control IMO selection via the flight computer keyboard. 
Thus, the crew would inhibit the redundancy management switching function during 
the period beginning just prior to taking the first measurements and continuing 
until after taking the second measurement. 

It is not necessary that the same IMO be employed for different rotation 
maneuvers. can be determined using a different MO. However, the appropriate 

alignment transformation must be enployed in the computation of There is 

a different T^® transformation corresponding to each IMO, unless the IMO stable 
irU 

platforms are aligned and torqued to the same inertial orientation. ' 
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EKROR ANALYSIS QF QSBITER EIGMyECTOR 


To begin th.e alignment transfer error analysis, the rotation eigenvector 
pointing error due to Orbiter UHJ errors ii?ill be determined. We recall that the 
eigenvector ^ is computed from rotation matrix C, as follows; 



Cfi 

*^12 


C = 

^^21 

^^22 

‘^23 


^^31 

*^32 

^^33 

c^ = 

^^^32 

" *^23’ 

^=13 


^ ~ ^*^13 ” ^^ 31 ^ ^*^21 



-h 


^ = kc 

The rotation matrix C is calculated as follows (derived in an earlier section) 


- vi - ’S' ’S 


,PC1 RL.t 
RL NB^ 


^ JC2 ^ ^ ^ PC2 ^ 

EL NB EL PCI EL PCI 


EL 

It is seen that is eliminated in the calculation of C. The errors associated 
with T^ are the Orbiter body-fixed IMO errors listed in Table 4 56 ^^ , 

^^COR^ ' geometrical misalignments of the Orbiter nav base and IMI. 

Thus, these errors do hot affect the accuracy of C. 


Introducing error perturbations yields: 




The errors 6T 


,PC2 

EL 


and are IMU readout random errors due to the giiibals and 


resolvers. They can be expressed as the effects of small error rotations of 








where 


page is 
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0 h ’ ° "Y2 

6Rj = Yl 5*1 Y2 ^ 

aj ^ _ j-^z ot2 

and ai„ 3 i, Ya “2» ^2» Y2 small angular (readout) errors about the X, 

Y, and Z coordinate axes respectively. 

Substituting yields: 

C + fiC = (I +6R2) (I + 

Since (I + 6X2)^ = I - 5 R^, we have 

G + 6C = (I + 6R2) C (I - 6R^) 
and, to first order, 

6C = 6R2C - C6Rj^ 

We will now choose the coordinate frame in which the eigenvector error 6^ 
will be ^alyzed. We choose the coordinate ffame such that the Z axis corresponds 
to the rotation axis of G; i.e., G represents a rotation around the Z axis through 
angle 0 . Figure 1 illustrates the situation. The normalization factor k and the 
eigenvector ^ are computed as derived above . 




Figure 1 . Eigenvector _d Error ^alysis Frame 


We now compute 6G in terms of the IMJ errors ai,6i,Yi. and 02, 62* Y2» Sub- 
stituting in the expression 6G = 6R2C - G6R^ yields 


50 




■ -0 

~Y 3 

2 2 


COS0 

-sin0 

0 ■ 


-Y sin0 
2 

<-Y COS0 , 
( 2 

3 

2 


Y 

2 

2 


sin0 

cps6 

0 


Y COS0 
2 

(-T sin0 
12' ‘ 

1 1 

-a 

2 


-3 

2 

1 

o 

CNl 

B 


0 

0 

1 


• < 
1 —3 cos0+a sin0 i3 sin9+a cos0 i 
[ 2 2 .2 2 ! 

0 


— 

— . 






1 

1 

. . 


cos 9 

- sln0 0 


0 

■^1 

3 

1 


^ 1 
-*Y sxn0 1 

L f 

“Y^cosSi gjCos0+ u^sin0 


C6R^ = 

sin0 

cos 6 ® 



0 

■^1 

— 

1 

Y COS0 . 

1 

1 

Y^sln0j 3isin0- aicos6 



0 

L. 

0 1 


-3 

1 

a 

1 

0 


i 

1 

a j 0 

1 ' 



Hence 



“ (y ~Y ) sin0 
2 1 

“(y -Y )cos0 , 3 

,21 I 2 

1 1 

- 3 COS0 - a sine 
1 1 

6C = 

(y “Y )cos0 
2 1 

I 1 

1 _(y -y )sin0 |-a 

'21 2 
< 

- 3 sine + a cos0 
1 1 


3 -3 COS0 + a sln0 
_ 12 2 

' 

I -a +3 sine + a cos0 } 

12 2 I 

0 


The error in jd, from the expression ^ = kc above, turns out to be; 

5d = k5c - (d • kSjc)^ 

Since 6_c = (50^2 - ^*^23* ^^13 ~ *^*^31’ ^*^21 ~ th^ It follows. 


6c^ = (a - a ) (1 + COS0) + (3 +3 )sin0 

^21 12 

6 c = (3 “3 ) (1 + COS0) - (a + a )sln6 

^21 12 

<Sc = 2 (y ~y )cos0 
^ 2 1 

Substituting 6c^ in the expression for 6^ yields: 

Sdjj = 


6dy = k6cy 


ORIGINAL PAGE IS 
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5d2 = 0 
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The result = 0 should not be surprising since jd by definition is a unit 
vector pointing in the Z axis direction. 

To con 5 >lete the error analysis of d, we determine the error statistics of 
6«L and 6dy in terms of the IMJ error statistics. Treating oi, gj, Yl and 
(X 2 , 32* uncorrelated random errors with zero mean yields the iollowing 

mean square errors: 

5c^ = (a^ + a^) (1 + cos0)^ + (3^ + 3^)sin^0 

X 2 a 12 

6c^ = (3^ + 3^) (1 + cos0)^ + (a^ + a^)sln^0 
X 2 J 11 

Also, 

6c^6c_ = (a,^ - + 3^ “ 3^) (1 +cos0)sin0 

X ^ 12 2 1 


From Table 4, the standard deviations of a , 3 


to 60, 


Thus, 


Y and a , 3 , Y are each equal 
1 2 2 2 


6c^ = 6cy = 4(60j^)^(l + COS0) 


6c^6cy = 0 

Computing 6d^ and 6d^ yields; 

6^ = 6d| == k26'^2 = 4(6^) ^(1 + cose) ^ 

4 sin^0 1 - cos 9 


Thus, the per-axis standard deviation of 6d is seen to be: 

-h 

dev6dj, = devddy = 60j^(l-cos0) 

We note the following; 

o For small 9, the per-axis error in ^ is very large. 

o When the rotation maneuver turns through 0 = 90 degrees, the per-axis 
error in ^ equals 50^, the IMD^ readout error. 

o When the rotation maneuver turns through 0 = 180 degrees, the per-axis 
error is 60 / iA2. This is the minimum error, as a fimction of 0. 

- 1 . * i . r ' ix 
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o ^ approaches 0 as 0 approaches 180 degrees. The particular calculations 
employed above to con^jute from the coiq>onents of C lead to a mathemati- 
cal singularity at 6 = 180 degrees . Such a conputation of c^ mechanize'd 
in a computer would "blow up" for 9 within a fraction of a degree of 180 
degrees due to conputer q^uantization, creating large errors in the compu- 
tation of We will see later that other equations can be employed to 
coopute from C. Euploying these different calculations does not change 
the error analysis results above. 
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ERROR MALYSIS 0? lUS EIGENVECTOR 


For the IDS, the measurement errors consist principally of the gyro mis- 
alignment, -scale factor and gyro drift errors- These errors are not random from 
measurement to measurement but are, rather, unknown parameters in their effects. 
The error expression derived for the Orbiter eigenvector does not apply to the 
IDS eigenvector. However, to first order, the IDS eigenvector error may be 
easily determined. 

IDS Gyro Scale Factor Error Effect 

The scale factor error causes an angular error proportional to the angular 
change sensed by the DID. The IDS eigenvector pointing error resulting from IDS 
gyro scale factor error is on the order of 6 (la) . This is shown immediately 
below. To siuplify the analysis, the eigenvector (axis of rotation) is assumed 
to lie in a plan containing the lDS*s pitch and yaw axes. We also assume non- 
redimdant IDS gyros, respectively oriented along the IDS pitch, yaw and roll axes. 


Let ky and represent the pitch and yaw gyro scale factor errors. 




Orbiter/IDS turns through 
angle 0 


Y Gyro 


where 

9 is the rotation angle of the alignment transfer maneuver 

(}i is the angle between the eigenvector and the pitch axis 

Then, the rotations 6^ and 0 sensed about the Y and Z axes are 

X Zt 

~ ky)6cos^ 

0^ = (1+ fc2)0sln(f> 


The eigenvector angular error 5^ is calculated as follows. 


tan(4> + 60) — 


(1 + k ) 0sin0 


(1 + ky) 0COS0 



This results (to first order) in 


50 



sin20 


3 ^ 



Notice that is Independent of 0 


Setting 




A 





^ y{<zUis 

siAZ<f 

'ff 


So, with = 45 ppm (lo) , then 

Dev ij/ == 6.6 sin 2 ^ sS 

We note the error is maximum when “ ^5 degrees. If ^ = 0 degrees or 90 
degrees, the error effect is zero. The scale factor error (in these cases) does 
not affect the sensed axis of rotation. We ^,jj^ll assume that the IDS gyro scale 
factor error effect on IDS aligp^nt is 6 sec (la) . 

IDS Gyro Drift 

The effect of IDS gyro drift depends on the elapsed time during the alignment 
process. If we assume 10 minutes is required to co 25 )lete the alignment maneuvers, 
then the per axis attitude error caused by gyro drift is about 6 sec (la) . 

IDS Gyro Misalignment 

The gyro misalignment directly affects the sensed eigenvector pointing 
direction. Thus, gyro misalignment introduces 10 sec (la) per-axxs error into 
the eigenvector per-axis error. 


IDS Eigenvector Error 

The per-axis IDS eigenvector angular error is .the^ESS of the three principal 
error sources discussed above. The ESS value is 13^ sec (lo). The IDS eigenvector 
error does not depend on the amount of Orbiter/IDS rotation, assinning some minimum 
amount of rotation necessary to eliminate the effects of small, random error 
sources . 

Different IDS IMP Models 

At the time this report was written, Boeing was considering two procurement 
sources for the IDS IMD. The error models are somewhat different. The differences, 
however, do not significantly change the analysis results. 
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ALIGMMENT TEANSFER EBROR 


The allgnioent transfer error is the sunt of the Orhiter alignment error, 
Orbiter ginibal errors, and lUS Dffl errors. The alignment transfer error is 
confuted as the ESS of the individual errors below. 


Table 7. Per-Axis lUS Alignment Transfer Error 





Orbiter 

lUS 

lUS 



Orbiter 

Transfer 

Eigen- 

Eigen- 

Alignment 

Case 

Alignment 

Error* 

Rotation 

Angle 

Vector 

Error 

Vector 

Error 

Error 



la 

3a 

■■ 

114 sec 

90 deg 

53 sec 

13 s^ 

126 s^ 

6.3 mlti 

i^HIH 

(Alignment 
per OFT 
Procedure) 

180 deg 

38 ^c 

13 sec 

121 s^ 

6.1 min 


48 sec 
(Modified 

90 deg 

53 s^ 

13 s^ 

73 s^ 

3.7 min 

B 

procedure, 

biases 

removed) 

80 deg 

38 sec 

13 sec 

63 

3.1 min 


*Does not include Orbiter IMU drift error effects. 


The table above does not account for the Orbiter* s gyro drift error 
(.035 deg/hr., Icr) accumulating between the times of Orbiter alignment and the 
Orbiter/IUS alignment transfer. Below, Tables 2 and 3 (repeated from the results 
summary) show the effects of Orbiter gyro drift. The effect of restricting star 
measurements (Orbiter alignment) to the 4x4 degree EOV is also shown. 


ORIGINAL PAGE IS 
OF POOR QliAL^f 


56 

























Table 2. Case A; lUS ALxgiHuent, Orbxter Sensor Body-Fixed Mxsalxgnment 


Errors Not Removed* 


Elapsed Txme * 

Alxgnment Transfer 
Rotatxon Magnxtude 

lUS Per-Axxs Alxqnment 
Accuracy (3cr) 

10 X 10 deg POV 

4x4 deg FOV 

0 mxn 

90 deg 

6.3 mxn 

6 . 0 nixb 


180 deg 

6.1 mxh 

5 . 8 -mxn 



^ ^ 


20 min 

90 deg 

6.7 min 

6.4 m^n 


180 deg 

6.5 mxn 

6.2 mxn 





40 mxn 

90 deg 

7 . 6 mxn 

7.3 mq.n 


180 deg 

7.4 mxn 

7.2 mxn 

60 mxn 

90 deg 

8.9 mxn 

8.7 l^fl 



180 deg 

8.8 mxn 

8.6 min 


* Note: 


o Elapsed txme xs perxod between start of Orbxter alxgniiient and end of 
Orbxter/IUS Alxgnment transfer, 
o Angle between Orbiter alxgnment stars assumed 90 degrees 
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Table 3. Case B:IUS Alignment, Orbxter Sensor Body-fixed 
Mxsalignment Errors Removed * 


Elapsed Txme* Alxgniaent Transfer 


lUS Per-Axis 


Rotatxon Magnitude 


Accuracy 



10 X 10 deg FOV 


0 

min 

90 

deg 



180 

deg 

20 

min 

90 

deg 



180 

deg 

40 

mxn 

90 

deg 



180 

deg 

60 

min 

90 

deg 



180 

deg 


3.7 ^In 

3.1 min 

4.3 raxh 
3.7 mxh 

5.6 mxn 

5 . 2 min 

7.3 m^ 
7.0 mxn 


Alignment 

(3ti) 


3.4 mxn 

2.8 iiixh 



4.0 min 
3.5 raxh 

5.4 rnih 

5.0 mxh 


7.2 :^h 
6. 9 mxh 


*Note: 


o Elapsed txme xs as defined for Table 2 

o Angle between Orbiter alignment stars assumed 90 degrees 
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PART 4 


DIPECT lUS ALIGNMENT PROCESS 
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DIRECT lUS ALIGNMENT VIA COMBINED PROCEDURES 


In the sunnnary of Results, the combined procedure was given for Orbiter 
alignment, body-fixed bias removal, and lUS alignment via rotation manewers. In 
essence, the Orbiter/IUS rotates 180 degrees about the LOS to each alignment star, 
simultaneously recording Orbiter star measurements and IMU attitude measurements 
(both Orbiter and IDS) before and after each rotation. The IDS alignment accuracy 
resulting is better than the Orbiter 's alignment accuracy and, in fact, is inde- 
pendent of the Orbiter' s alignment. 

Rationale for the Combined Maneuver 


Imagine that the flight crew were able to orient the Orbiter so that the star 
image always fell exactly in the same location in the star tracker's FOV (e.g., 
directly on the tracker's boresight axis). Then, the Orbiter' s eigenvector 
produced by the 180 degree maneuver around the star LOS would be exactly colinear 
with the star LOS. In the IDS computer, the IDS eigenvector, although expressed in 
the unknown IDS inertial coordinate frame, would also point exactly at the selected 
star. 

Since the star's inertial coordinates are known in the desired inertial 
navigation coordinate frame (e.g., M50) from star catalogue data, it becomes a 
single matter to calculate directly the IDS alignment transformation relating the 
two coordinate frames. This would be done without regard to the Orbiter' s IMD. 

The only errors bearing on the IDS alignment achieved thusly would be the Orbiter 
star tracker measurement error and the IDS IMD errors. The resulting IDS alignment 
error would be about 1 min (3o) per axis, which is considerably less than the 
Orbiter' s alignment error. (This will be shown later.) 

In the actual situation, we do not expect the star image to always fall 
exactly in the same FOV location. However, the image will always fall somewhere 
in the FOV and hence the eigenvector produced by the 180 degree rotation maneuver 
will not deviate by more than 5 degrees maximum per axis- from the^tar LOS unit 
vector. 


Now, it would seem since the measured star LOS unit vector and the rotation eigenvector 
are almost colinear (within a few degrees) that information is available to essen- 
tially eliminate Orbiter IMU errors from the IDS IMD alignment problem. It turns 
out that this is the case, as discussed next. 


Computational Approach Taken 


We begin with the fact that the star LOS unit vector and the rotation eigenvector 
point essentially in the same direction. The following steps are then taken: 


( 1 ) 


Compute the difference vector (small) between the star LOS unit vector and the 
eigenvector, in Orbiter stable numbers coordinate. 


Measured - 
Star LOS 
, Unit vector 
‘ (Orbiter Star 
‘ Tracker) 



Computed difference 


Measured Orbiter 
Eigenvector 
(Orbiter IMD) 


Orbiter Stable Member Coordinate Frame (Present Cluster Frame) 



(2) Transform the difference to the lUS unknown frame using transformation' Tp^ 

(present cluster to unknown frame) . This transformation is determined from 
the Orbiter and lUS eigenvectors expressed in Orbiter IM0 stable member and 
IbS IMD unknown frame coordinates respectively. , 

(3) Add the transformed difference to the lUS eigenvector, to con 5 >ute the star 
LOS unit vector in the lUS unknown frame coordinates. 

Transformed difference 

measured XUS 
eigenvector 
(XUS XMD measurements) 


(4) Use star LOS unit vector (U frame) and corresponding star catalogue data 

(e.g., M50) to compute the desired XUS alignment transformation (e.g., T^ ). 

This completes the alignment determination. 

We see that the Orbiter IlfiJ errors enter only via the transformed difference 
vector. Because this difference is small, less than a tenth of the magnitude of 
the unit vectors, the effects of Orbiter IMU measurement errors for all practical 
purposes are eliminated. This will be shown in the error analysis. Also, it is 
seen that Orbiter XMU alignment or lack thereof is not a consideration. All Orbiter 
measurements are expressed in stable member coordinates; knowledge of the alignment 
of the IMU stable member relative to the M50 frame is not required. 

XUS Alignment Calculations , Combined Procedures 

Let us now proceed through the steps listed above, enploying the needed calcu- 
lations. We will assume, without any loss of generality, that the desired XUS basic 
navigation frame is M50. 

'i 

We will express Orbiter and XUS inertial attitudes as 3 x 3 matrices of body 
axis cosines. The columns of each matrix are formed by vehicle X, Y, and Z axis 
unit vectors respectively. For the Orbiter, the X, Y, and 2 axes (unit vectors) 
are expressed in the XMU stable member frame. This matrix exists in the Orbiter 
flight computer and is denoted T^g, the transformation from the nav base reference 
axes to the present cluster (stable member) frame. Xn the XUS flight conputer 
a similar quantity would exist (probably expressed as quaternions) relative, to the 
XUS unknown inertial coordinate frame. For this analysis we will denote the XUS 
inertial attitude as the 3x3 matrix U. 


Computed star 
LOS unit vector 


XUS Unknown Coordinate Frame 
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Coiq>ute 3 X 3 E matrix from star 1 and star 2 M50 coordinates (Star catalogue 


data) . 


I . , 

A e^~ unif ( e, X 6;,) 




e. 




02 (Mfo) 


r 1 /I 


5'^2S.r 2 


o Sight on the first star and take the following measurements 

star LOS unit vector, in stable member coordinates, first 

measurement first star. 

Orbiter attitude, first measurement first star, stable member 

coordinates . 

TJ__ lUS attitude, first measurement first star, lUS unknown 

frame. 

o Rotate around the first star LOS 180 degrees and repeat the measurements 

Star LOS unit vector, in stable member coordinates, second 

measurement, first star. 

PC12 

Orbiter attitude, second measurement first star, stable 

member coordinates 

lUS attitude, second measurement first star, lUS linknown 

frame. 


o Conq>ute with first star measurements 
^ = unit + %2^3 

<=1 = 

- ®12 


(bias removal) 



o From Ct and F-, , coiq)ute (gee Appendix A) first rotation eigenvectors 
ll and 

^ Orbiter rotation eigenvector, first 180 degree rotation, 
stable member coordinates 

g-j lUS rotation eigenvector, same rotation as for lUS 

taiknown frame 


Averaging and m^^ removes Orbiter body fixed biases from the star 
measurements (Appendix D)i is the orthogonal transformation, in stable 
member coordinates, linking the initial and final Orbiter attitudes 
relative to the 180 degree rotation around the star LOS. 

is the orthogonal transformation, in lUS unknown frame coordinates 
linking the initial and final lUS attitudes relative to the same 180 
degree rotation. It is from and F, that the Orbiter and lUS eigen- 
vectors and ^ are determined, as ^own in Appendix A. Although 
expressed in different coordinate systems, ^ and point in the same 
direction in inertial space, given that the lUS position relative to 
the Orbiter is unchanged from measurement to measurement. 

t 

o Sight on the second star and take the following measurements 
(simultaneously) ; 


m 9 i star LOS unit vector, in stable member coordinates, first 

measurement second star 


tJC21 

•^NB 


Orbiter attitude, first measurement second star, stable 

member coordinates 


attitude, first measurement second star, ITJS unknown 

frame 

o Rotate around the second star LOS 180 degrees and repeat the measurements 


^2 “■ 


star LOS unit vector, in stable member coordinates, second 

measurement, second star 


t^C22 _ 


Orbiter attitude, second measurement second star, 
member coordinates 


stable 


^22 attitude, second measurement second star ITJS unknoxim 

frame 


o Compute with second star measurements 
22 = mit |%(m2]^ + ^ 2 ^ 

C2 - 

22 [?2l] ^ 


= U, 
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o From and F 2 » compute second rotation eigenvectors ^ and jg 2 C&ppendix A) . 


-S2 


Orbiter rotation eigenvectiouj second 180 degree rotation 
stable member coordinates 


— lUS rotation eigenvector, same rotation as for ItlS unknown 
frame 


o Form 3x3 D matrix from Ai and ^ 
^2 ~ unit X 


d2 = ^2 ^ 


H - [ii |i2 iis] 

o Form 3x3 G matrix from and ^2 - 

_g2 = ^it (^ 3 ^ X 

^2 "" ^ ^ £1 
G = 13.2 1 ^ 3 ] 

o Confute difference vectors and A^ 2 * 


Ail “ Si " ii 


Ai 2 = S2 - S2 

The A_d^ and Ad? are the differences between the star LOS unit vectors and 
the eigenvectors, in stable member coordinates. 

o CoD^jute transformation TpQ relating the Orbiter IMO stable member (present 
cluster) frame to the ITJS unknown frame. We know by definition that 
transforms d^^ to ^2 £, 2 » is Thus TpQ transforms D to G. 

G - T«c B 

Solving for yields 

= GDt 

o Transform differences A^j^ and A ^2 the ITJS Tjnknown frame (call the 
results A^j^ and A^ ) 

Ai 


- d^ 


142 - '^■^2 


hh 



o Confute star LOS 'unit vectors and li^ 

= mit 

= unit (A £^2 + ^ 2 } 

The above yields conputed star 1 and star 2 LOS unit vectors in lUS 
unknown frame coordinates 

o Form 3 x 3 H matrix from II 2 * 

= unit X h^) 

h* = h 3 xh^_ 



Compute the lUS alignment transformation 


Compute the lUS alignment transiormatlon I'lj’ relating the IDS unknown 
frame to the M50 coordinate system. We know that by definition 

transforms h^^ to to ^3 to ^ 3 . Thus, transforms E to E. 


E = T^^H 

D 

Solving for 3 Tields 


t“50 ^ 


T^ is the desired result of the alignment procedures. The actual 
alignment is carried out sTibsequently by transforming IDS inertial atti- 
tude D from IDS unknown frame coordinates to M50 coordinates iising T^®. 
This is a one-time computation. The IDS strap down inertial system will 
then proceed to update IDS inertial attitude in the M50 system. 
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PJOIT 5 


ERROR MALYSIS 


DIRECT IDS ALIGNMENT 
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EKROR MALYSIS OF DIKECT lUS -ALIOfflENT 


From the equations above, we see that the lUS alignment error depends on 
the error in h, the confuted star LOS unit vector e^ressed in lUS miknown frame 
coordinates. The equations are: 

h = As. + ^ 

Ai = 

Ad = m - d 


where 

m star LOS unit vector, average of two measurements. Orb iter stable 

member coordinates 

_d, ^ — Orb iter and lUS jointly sensed eigenvectors, in Ofbiter stable 
member and lUS unknown frames, respectively 

Transformation from Orbiter stable member frame to lUS unknown frame 

We will assume the star catalogue data to be error free. We will also assume, for 
the present, that the two alignment stars are 90 degrees apart. 

The error equation for li is : 

6h = 5^ + [6T^J Ad + 

Determining the errors fi^,6T^^,and 
error. 

Error fig 

The error 6g^ is due entirely to the IDS IMD. This error has already been 
determined to be 13 s^ (15) per axis. 

Error 


The per axis error STp^ has already been determined. It is the ESS of the 
Orbiter eigenvector error (38 lo) and the IDS eigenvector error (13 ^^,la), 

which is 40 sec. The effect of 5 T^q is to rotate 40 sec per axis. 

* u 

Si nce A^ is small' con^ared to ^ (a^ tenth or less) , the error STp^A^ is less 
than 4 sec (la) per axis. The smaller M, the smaller the error effect. 


6A^ will allow determination of the IDS alignment 
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Error SAd 


Both m and jd incorporate the same IMtJ measurement. In addition, m incorporates 
the star tracker measurements. It would appear, since the IMD measurement errors 
are common to both m and jd, that in forming the difference ^ the IMCF errors might 
cancel, leaving only the star tracker error effects. This is indeed the case, as 
demonstrated below. 

We will build on the error analysis already cou^jleted determining the error 
in the eigenvector It_was determined that 

== Ain{^®2 "V ®] 

Sd^ = 0 
where 

^ 1*^1 ^ ^ axis angular measurement errors due to Orbiter IMD gimbals 

and resolver, first measurement (prxor to 180 degree rotation) 

”2*^2 ^ ^ SCS 3 .S angular measurement errors, second measurement (taken 

after the 180 degree rotation) . 

The above results were determined in a coordinate frame in which the Z axis 
is oriented along eigenvector Employing the measurement errors 0j, & 2 * 

the averaged star tracker measurement error 6m will he calculated. This will be^ 
followed by the calculation of 5 Ad. To simplify the analysis, we will assume that 
m also points in the Z axis direction. To first order, this will not affect the 
error analysis since m and ^ are only a few degrees apart. 

m is the average of two measurements (in stable member coordinates) which 
will be labeled here as m^ and 32 • Thus 

m = h (^ + 

6m = % (6^ + 6^2) 

Averaging m^ and m, removes the Orbiter body fixed sensor biases from m (as 
explained in Appendix D) . We need consider only the random errors. These are 
the star tracker measurement errors (denoted 68gy in the error model. Table 4 ) and 
the IMD gimbal/ resolver error (denoted Table 4 ). Thus, the errors 6mj^ and 

6^ are the random angular error effects of the Orbiter IMD and star tracker 
measurements . Dsing nomenclature previously employed (page 30 ) , the errors can 


be written: 


= 6R^m^ 

ORIGINAL PAGE IS 
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where and 6R2 .are now the rotation error matrices (3 x 3 ) due to the Orbiter 
and star tracker errors, first and second measurements respectively. 
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The per-axls measure m ent errors corresponding to 5R^ and are written below. 
First measurement; 

X-axis angular error = 

Y-axis angular error = 

(z-axis angular error we ignore since m^ = (0, 0, 1)) 

where 


HHJ errors (radians) 

^3’ ^3 star tracker errors (radians) 

Second measurement: 


X-axis angular error = 02 + 0 ^ 

Y-axis angular error ~ ^2 

(z-axis angular error we ignore since ~ 0, l)j 

where 


« 2 » ^2 IMB errors (radians) 

u^, B^ star tracker errors (radians) 


Thus 


6R^ 


0 

0 


-(3,H-B3) 


0 

0 


(Ui+Ub) 


(3^+ B 3 ) 

-(a^+ 03 ) 


5R. 


0 

0 




0 

0 


Ca2+a^) 


(32^- 3^) 

■(U2+ «4) 


ORIG^AL PAGE IS 
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We have, to first order. 


6m = + 6m2) = %(6R^ + 6^2)2 

Thus 


(3^+62+33+3^0 

6 m ~ \ , 

~ -(a^+a2+a2+a^) 

^ 0 - 
Now calculate 6 Ad = 6 m * 5^. 

5 Ad^ = 35(3^+62+33+3^) - 2 sL e 

fiAdy = -% ( 0 j+a 2 +a 3 +a^) - ^ ^ -3^) (1 + cos 0 ) - (a^^ + 02 ) sin ^ 

AAd^ = 0 


Combining terms yields 


6 

6 Ad^ = - 

6Ady = - 



”3 -^“4 
2 


(“2 ” 02 ) (1 + cos 9) 

2 sin 0 

(3j^ - 62 ) (1 + cos 9) 
^ 2 sin 0 


6Ad^ = 0 


The mean square errors are; 



where it is assumed that ,3, ,« 2 » 3 p, 03 , 60 ,( 2 , , 3 . are uncorrelated with zero means. 
From the error model. Table 47 we have ^ 



Thus 



dev(6Ad^) = dev(6Ady) = 





dev(6Ad2) = 0 
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The eigenvector rotation angle is 180 degrees. For 6=180, (1+ cos 0>/ 
(1-C03 8) =0; the per axis error in $Ad, to first order, is due only to the star 
tracker error The IMO errors ''cancel out". 


dev(SAd^) = dev(5Adj) = 


d8 


ST 


vr 


(0 = 180 deg) 


The error model (Table 4) gives 66g the value 42 (1 a), ffe will also consider 

the situation where star roeasureioenes are restricted to the central 4x4 degree 
FOF of the star tracker, which it is understood halrves the error. Thns: 


Table 8. Ad per— Axis Error Value 


FOV 

1 

8 = 180 deg 

10 X 10 deg 
4 X 4 deg 

30 ^ (15) 
15 s^ (15) 


It iS recalled that 6Ad effects lUS alignment accuracy in the tern (5Ad). 
Since Tp^ is an orthogonal transformation, the per-axis error tS^ (5nS) remains 
the same as given in Table 8 immediately above for 6Ad. ~ 
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lUS jJLIQBffiNT ERROR, DIRECT ALIGNMENT 


From sibove, it was determined that the error in h, the coiqmted star LOS 
unit \fector expressed in IDS unknown frame coordinates, is 

6h = 6^ + [sTp J Ad- + Tpj,(6Ad) 


A rigorous treatment at this point would account for the correlations between 
(l)6Ad and StJc and (2) 5 T?c and %, since ^ and ^ were en^loyed to calculate 
^?C* purposes, we ignore the correlation because the major error sources 

and 6A^) are uncorrelated. These errors are dxie to the star tracker and lUS 
DflJ respectively. 

We consider the selected alignment stars to be 90 degrees apart. Hence, 
the per-axis IKS alignment error equals the per-axis error in SIi (Appendix C) . 

We thus calculate the IDS alignment error as the ESS of the conqjonent errors. 

The results are presented immediately below. 


Table 9. lUS Alignment Error (Alignment Stars 90 Degrees Apart) 


Star Tracker 
FOV 


.4 .a 

Tp^.(SAd) 

IDS Alignment Error 

lO 

3<J 

10x10 deg 

13 s^ 

4 sec 

30 s^ 

33 ^c 

1.6 min 

4x 4 deg 

13 sec 

2 sec 

15 

20 s^ 

1.0 min 


lUS 

IMD 

Error 

Effect 

Transfor- 

mation 

Error 

Effect 

Star 

Tracker 

Error 

Effect 




If the alignment stars are not 90 degrees -apart, then the alignment degrades 
somewhat, as explained in Appendix C. For alignment star angles within 90 + 30 degrees, 
the average degradation is less than ten percent. 
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PAST 6 


STRAWMAN MECHANIZATION 


DIRECT IDS ALIGNMENT PROCEDURE 
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DMA, BirEKFACE 


■*r 


The equations developed above were exaniined for on-board implementation in 
case such were considered. A brief investigation sRowed implementation should be 
straightforward. The description below represents a possible implementation 
approach. Figure 2 presents a block diagram of the Orbiter/IDS data processing 
interface . 

Software Mechanisation 


It is assumed that the alignment equations would be processed in t’hp lUS 
computer. Orbiter data would flow from the Orbiter GEN&C flight computer to the 
Orbiter Systems Management (SM) flight computer to the lUS flight computer. 

A small software program would be required in the SM computer to control 
the alignment transfer. This program would accept keyboard inputs from the crew, 
notify the IDS computer of the iuqjending alignment, monitor the GN&C computer for 
the start of each data collection period, signal the IDS con^uter to take IDS 
data at the appropriate times, receive and transfer Orbiter data to the IDS flight 
computer, and provide CRT displays for crew control. 

The only new software required in the GN&C computer would be a flag in the 
COMPOOL data base, set and reset at the initiation and completion of each Orbiter 
data collection period (lasting 3.2 seconds). The SM computer would monitor this 
flag every 160ms, when in the IDS alignment mode. 

Operational Sequence 

The IDS direct alignment operational sequence would be the following; 

o Crew commands Orbiter IMD in-orbit alignment via the GM&C computer. 

o Before the first star sighting is taken, the SM computer is placed in 
the IDS alignment mode by keyboard command. 

o Thereafter, four collections of Orbiter and IDS data sets take place, 
as the Orbiter maneuvers and takes four star tracker sightings (on two 
stars) in accord with the direct IDS alignment procedures of this docu- 
ment. Data would be automatically taken and transferred to the IDS flight 
computer. 

o The SM computer would provide appropriate outputs to the CRT display for 
crew monitoring of the alignment process. 

Data Interfaces 


Orbiter data supplied to the IDS computer would be that indicated below. 
This data is located in the Orbiter GN&C flight coDq>uter COMPOOL data base. The 
nomenclature is that defined in the IMU SOP, reference 4. 


54 




Figure 6. ON-ORBIT DATA COLLECTION HARDWARE INTERFACES 















o I_STAI^SEL1, I_STA^SEL2. . . .unit LOS vectors from the center of the earth 
to alignment stars #1 and #2 respectively, M50 coordinates. 

o I_STAE|^SELC(I,J) Orhiter star measurement LOS mit vector, in present 

cluster (stable menibers) coordinates, IMO (1, 2, or 3), for the 

star (1 or 2) . This quantity is the averaged result of 21 star 
tracker measurements taken 160 milliseconds apart, representing 
data spread over 3.2 seconds. 

PCf J") 

' 3x3 transformation matrix, Orhiter nav base reference axes 

to present cluster coordinates. This quantity is the averaged result of 
21 Orbiter IMCJ readings, HHJ (1, 2, or 3), taken at the same time as 
the 21 star tracker measurements above. 

o Time tag... data time tag for measured quantities above. 


GN&G Computer 

The GN&C computer would set a flag in COMPOOL when the star tracker measure- 
ments begin. Specifically, the flag is set when star tracker software DAT^FILR 
routine is called by the GN&C operating system (FCOS) . 


Orbiter SM Computer 

The ITJS ali g nm ent process would be controlled via the Orbiter systems manage- 
ment (SM) flight computer. The SM computer would alert the lUS computer that the 
alignment transfer is about to take place. SM keyboard entry would determine 
which Orbiter IMU (J=l,2, or 3) would provide data to the lUS. 


The SM con 5 >uter has the capability of accessing SM selected GN&C COMPOOL data 
via the ICC_SSIP software module. The GN&C data is transferred to the SM computer 
via inter-conqjuter channels (ICG) . This data is then quickly transferred by the 
SM con 5 >uter to the lUS computer via payload Interface. 


The SM computer would do the following basic things: 
o Cyclically (every l60 ms via ICC) receive DATA_FILR flag, 
o When flag is set, command the lUS con 5 >uter to take ITJS attitude data, 
o Request, receive, and transmit alignment data from GN&C to ITJS. 


( 1 ) 

( 2 ) 

(3) 

(4) 


I_STAR_SEL(I) 

I_STAR_SELC(I,J) 

T,PC(J) 

NB 

TIME TAG, for (2) and (3) 
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IPS Coniputer 

The lUS coii 5 >uter woiild collect the Orhiter/XUS data and perform the align- 
ment calculations, doing the following basic things: 

o Initiate IPS body attitude data collection (body attitudes and time 
tags) following receipt of command from the SM conqmter. This data 
would be collected, say, over a two second period every .2 second for 
filtering and interpolation purposes. 

o Receive Orbiter alignment data from the SM computer. 

o Inteirpolate the IPS data collected to the Orbiter time tag point. 

o Calculate the IPS IMP alignment matrix. If the IPS basic navigation 
coordinate frame is other than M50, the IPS coiqjuter would need a fixed 
3x3 transformation matrix to transfer I_STAR_SEL (I) to the desired 
coordinate frame. 

Timing Considerations 

The mechanization approach above has no critical timng requirements . Data 
staleness is not a problem since I_STA^SELC(I,J) and are time tagged. 

I__STAI^SEL(I) are constants (vectors) and they form part of the GM&C computer 
I-Load (pre-mission data load) . 

The approach above does assume that the Orbiter and IPS flight computers 
operate with a common time base, since time tags are involved. It is understood, 
tmofficially, that there will be a hardwire connection between Orbiter and IPS 
master timing units such that the respective time bases will be significantly less 
than a millisecond apart. This accuracy is entirely adequate. 

The only real-time requirement is that each IPS data collection period 
(2 seconds suggested above) fall within the Orbiter 3.2 second data collection 
time span. Thus, the IPS data collection should begin no later than .75 seconds 
after the Orbiter data collection begins. The brief investigation conducted to 
date indicates that the IPS coiqjuter can be notified about . 33 seconds (maximum 
delay) after the Orbiter data collection process is Initiated. Thus, there 
appears to be plenty of time for the IPS cor 5 >uter to Initiate its data collection 
process. Subsequent processing of the Orbiter and IPS collected data by the IPS 
computer is non-cyclic and non- time critical. 
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PART 7 


CONCLUSIONS 
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CONCLUglON 


The basic conclusions of this report are: 

o lUS alignment can be achieved with an error much less than 6.3 

(3a) , independent of the Orbiter IMD alignment error. This is accom- 
plished by combining the OFT in-orbit IMtJ alignment procedure, modified 
to remove Orbiter sensor biases, and the lUS alignment transfer pro- 
cedure. The accuracy achieved thereby is estimated to be about 1.6 
m^ (3a) . 

o Accuracy can be in^roved to about 1.0 min (3a) by restricting star 
images to the central 4x4 degree star tracker field of view. 

o The recommended rotation maneuvers do not appear difficult to perform. 

o On-board in^tlementatlon of the recommended lUS alignment approach^ were 
such to be considered, appears to be straightforward. Impact to the 
Orbiter/IUS interface appears to be minimum. 
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APPENDIX A 


ORTHOGONAL TRANSFORMATION 
EIGENVECTOR COMPUTATION 


A-1 






APPENDIX A 


Determination of Orthogonal Transformation Eigenvector 

Let C be a 3 X 3 orthogonal transformation (rotation matrix) . Let x be an 
eigenvector of C. Since the eigenvalue of an orthogonal transformation is 1, 
then by definitxon 

Cx = X 


or, 

(C-I)x = 0 


Writing the equatxon above in component form yields 


‘^ ir ^ ^12 ^13 


^1 


0 

^21 ^23 


^2 

= 

0 

C 31 c ^2 


x3 


0 


The above is a Ixnear system of 3 homogeneous equatxons in 3 unknowns (x^^, x^, x^) 
of rank 2. The solution comes about immediately from a well-known property 
of determxnants . 


Solution 


Let ]be an nxn matrix with elements a^j . Let the cofactors of a- • be 
denoted A^j^ A property of determinants (termed the Laplace development is 





Denoting the matrix of 


cofactors as 



we can write the above as 




where I is the nxn identify matrix. 


det 



I 


det 


Applyxng the above, we letja^J - C - I. 
a^jj= 0. Thus, xn this particular case. 


Sxnce C - I is of rank 2, then 




= 0 


^ QVAurr 


A-2 



It is seen that any column of (i.e., any row of A^j ) is a solution to the 

linear homogeneous equation x = 0. Thus we can write the solutions 


= kA._ 

il 


X 


2 



X 


3 



where k is any arbitrary constant, and i = 1, 2, 3. 

The theory of linear homogeneous e<luations tells us that when the rank = n - 1, 
there exists exactly one independent solution. Therefore, the ratios xj/x^ are 
uniquely deteirmined; the solutions obtained for different values of i above are 
identical in terms of the ratios, although the computations are different. 

We now write down the three alternative ways to compute the solution of 
(G - I)x - 0, corresponding to i = 1, 2, 3. Each alternative consists of com- 
puting the appropriate cofactors . The cofactor of an element of C-I is determined 
by striking out the row and colu3m containing the element, forming the determinant 
of the 2x2 submatrix remaining, and assigning the polarity (-1)^ Let k = 1. 

Alternative 1 (let i = 1) : 

^1 ~ ^*^22 “ ^^^*^33 “ ^23 ^^32 

^2 " " ^^21 *‘'^33 “ ^23 ^31 

^3 " ‘^21^^32 “ ^^^22 " 

Alternative 2 (let i = 2) : 

^1 " ^12 ^*^33 " ‘^ 13'^32 

^2 = ^‘^ll " ^^^33 " “ '^13^31 

^3 " "‘^32^‘^U " + *^12 '^Sl 

Alternative 3 (let i = 3) : 

^1 " ^12*^23 “ ‘^13^‘^22 “ 

^2 ^ "‘^23^^11 '^13^21 

^3 " <^11 " ^^=22 - - ^12^^21 

The solution's x above must be normalized in order to form eigenvector ^ of 

unit^lerigthjy ,,Tfie normalization factor k = [x | ^ is, in general, different depending 

on th'e tcomputatidnal alternative chosen, 

V If. ‘ 



Of particular interest are the possible singularities which might arise 
during the computations. In this case, the matrix G given is such that under 
the computational alternative chosen, the trivial solution x " 0 results. This 
occurrence was noted earlier in the text of this report, page 33. However, since 
exactly one independent non-trivial solution exists, at least one of the compu- 
tational alternatives above must provide finite, non-trivial x. 


Take the situation from page 30 where C represents a rotation 
coordinate axis. 


C = 


COS0 -sin9 
sin0 COS0 
0 0 


0 

0 

1 


aroxmd the Z 


When 0 = 180 degrees, the conputational method employed (which is different from 
the three alternatives above) yielded x = 0, which is a solution, but not the 
non-trivial solution sought. Let us apply the three alternative computations of 
this section. The results are as follows: 


Alternative 1: x “ ^ (trivial solution) 

Alternative 2: x = 0 (trivial solution) 

Alternative 3: x = (0, 0, 4) 

Alternative 3 yields unit eigenvector^ = (0, 0, 1), which is the proper solution 
since we postulated C to be a rotation about the Z coordinate axis. 

Computational Procedure 

A sure way to avoid the singular solution 2 c = 0 is to coiq)ute x^ using each of 
the three alternatives. At least one solution will be non-trivial. If at least 
two solutions are non-trivial we could choose the solution x having the largest 
magnitude |x[. This minimizes the effects of computer roundoff errors. 

Which Way Does x Point? 

The non-trivial solution -x satisfies (C - I)x = 0 as well as the non-trivial 
solution X. Both x and -x lie along the axis of rotation (eigenaxis) represented 
by G but point in opposite directions. In applying the computational procedure 
discussed above, the question arises "which way along the eigenaxis does the 
computed solution x point‘d" The answer is inq)ortant, because we want the jointly 
sensed and computed eigenvectors (Orbiter and lUS) to point in the same direction 
in xnertlal space, not opposite directions. Also, when combining the Orbiter and 
lUS alignment procedures, we want the jointly sensed eigenvectors to point toward 
the alignment star, not away from it. 

Since the solution to (C - I)x = 0 is not uniquely determined as to its 
direction along the eigenaxis, an additional step must be taken. A relatively 
simple procedure resolving the question of "which direction?" i& described below. 
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Tot the sake of simplicity, it is assumed that the lUS is mounted in the Orbiter 
payload bay with its body axes (roll, pitch, yaw) roughly parallel to the Orbxter’s 
body axes (say within 10 degrees). The' basis of the procedure is to (a) transform the 
jointly sensed eigenvectors ^ (Orbiter stable member frame) and ^ (lUS unknown 
frame) to the respective body axis frames and then (b) make some con^tarxsons . Since 
the Orbiter and lUS body axes are assumed roughly aligned, the transformed eigen- 
vectors, in respective body axis coordinates, should be roughly eq^ll valent. The 
magnitude of the difference between the transformed eigenvectors should be signifi- 
cantly less than 1. If not, a reversal in direction is required. 

The following quantities are involved, all of which have already been employed 
in the main text analysis of the lUS alignment. 

^ _ Orbiter eigenvector, stable member coordinates 

^ XUS eigenvector (jointly sensed with ^) , lUS unknoxra frame 

coordinates 

FC2 

Tjjg transformation, Orbiter nav base to stable meniber coordinates, 

at time of 2nd measurement taken at completion of eigenvector 
rotation maneuver — the coltanns of this matrix are the Orbiter 
body axis direction cosines relative to the stable member frame. 

U 2 - - - lUS body axis matrix at same time point the columns of this 

matrix are the lUS body axis .direction cosines relative to the lUS 
unknown frame 


m - - “ - unit vector pointing toward alignment star, Orbiter stable member 
frame 

Ad - -- - A^=m-^; hence I A^[ « 1 


Eigenvector Direction Resolution Subprocedure 


The subprocedure begins at the point where the difference ^ = m - ^ is confuted. 

(1) Correct the direction of if required. 

If: Ij set .d 

Otherwise: ^ already has correct direction 

(2) Transform jd to Orbiter nav base coordinates 

i’ i 

(3) Transform ^ to lUS body axis coordinates 

(4) Correct the direction of if required. 

If: |d_’ - ^\> 1, set 

Otherwise: £ already has correct direction 
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This completes the subprocedure; ^ and ^ now have correct orientations in 
therr respective xnertial frames. If the star measurement is not xnvolved 
(separate Orbiter alignment and Orbiter/XTJS alignment transfer ) , simply start the 
subprocedure at step (2) . 

If it turns out that the Orbiter body axes and lUS body axes in the payload 
bay are not roughly parallel, the transformation relating these axes can be 
easily determined (prior to flxght) within the accuracy needed (say, within 10 
degree accuracy). This transfonnation, call it tJ]^, would be used to transform 
jd' to where d" is the Orbiter eigenvector expressed (roughly) in lUS body axxs 
coordinates. The instead of , would then be used xn step (4) above. 



APPENDIX B 


ORBITER IMO AND STAR 


TRACKER ERROR MODEL DISCUSSION 





ERROR MODEL 


Each of the three IMti’s consists of a stable platform interconnecred to 
its case through four gimbals . Each IMU case is mounted on a common baseplate 
called the navigation base. Also mounted on the navigation base are two star 
trackers. The body fixed orientation of each sensor case relative to the navigation 
base is carefully measured. 

The principal function of the stable platform is to maintain a fixed, 
known orientation in inertial space, regardless of attitude changes by the Orbiter 
vehicle. Gyroscopes on the platform sense angular motion of the stable platform 
and act through servo loops and giaibal motors to maintain the platform orientation 
fixed in inertial space. The orientation of the Orbiter’ s navigation base relative 
to the stable platform is determined by simultaneously measuring the four IMtl 
gimbal angles (gimbal resolves). These angles, plus the known orientation of the 
IMU case relative to the navigation base, allow the determination of the Orbiter ’s 
attitude relative to inertial space. 

Complete error models of the IMU and star tracker address eveiry error 
source important to overall operational capability. An example of a complete, 
detailed error model is Reference 25. For the purpose of analyzing in-orbit alignmeni 
and alignment stability, only a subset of these error sources is needed. These 
error sources turn out to be the mcompensated errors associated with (1) sensor 
geometrical orientation errors relative to the navigation base, (2) IMU gimbal 
and resolver errors, (3) gyro drift rate errors, and (4) star tracker measurement 
errors . 


Table B1 presents the error model employed in this analysis. The primary 
data source is the error budget presented by Rockwell International (RI) at the 
15 March splinter meeting. These errors have been further broken down (Table B2) 
into error types in accord with reference 5. Discussion of these errors follows. 
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TABUE B1 


ALIGNMENT ERBOE MODEL 

ERROR SOURCE SYMBOL ^ YALUE/AXIS (la) 


Orbit er IMU 


Gyro Drift 


Nav Base Ref. to Mounting 
Pads 

6e 

'^IPN 

IMU Case to Pads 

®^CMP 

Case to Outer Roll Gimbal 

^®C0R 

Non-Orthogonality | 
Resolvei J 



Star Tracker 

Horizontal, Vertical *^®ST 

Measurements 

Tracker to Nav Base Ref. <Se.„rn 

INdJ. 


.035 deg/hr 

60 

20 

28 


Body 

Fixed 

Biases 


30 sec j 53 sec 
^c J RSS 


1 68 sec 
RSS 

sec 


60 sec I Body Fixed Bias 


RSS (not including gyro drift) e 


llU s^ (1 c) 
5.7 3^ (3 o) 


* Note. These symhols are employed in the error analysis. 


B-3 



TABLE B 2 


BREAKDOWN OF IMU CASE TO STABLE MEMBER ERRORS 


ERROR SOURCE 


ERROR VALUE/AXIS (la) 


Pxtch. to Outer Roll Gimbal 

Won-Orthogonality 

c Resolvers — lj.l{ 


Component ■ 
Errors 


Offset 


Random bias 


Sinusoids 


1 speed 

7.6 

8 speed 

19 

9 speed 

k .2 

16 speed 

20 


Quantization (20 s^) 


30 sec 
12 s^ 

29 

6 sec 


ESS subtotals sec 


c Case to Outer Roll Gimhal 

Geometrical misalignment 


ESS (I0) 


28 sec 


60 sec 




b-4 



The Geometrical Orientatxon Errors 


The body fixed geometrical misalignments of the IMU's and the star trackers 
adversely affect the precision with which the star sightxng vectors can be 
transformed xnto the IMU stable platform inertial reference frame. 

These errors stem from (1) the precxsxon with which the sensor case orienta- 
tions on the navxgation base can be optically measured prxor to lift-off and (2) on 
flexing of the navigation base in orbit. The flex xs due principally to the heat 
load imposed on the navigation base by the IMJ’s and star tracker's power supplies. 
During the period of time between the Orbiter star alignment and the alignment 
transfer to the ITJS, the flex can be considered constant. 

An addxtional body fixed misalignment (per reference 2) is the IMU case-to- 
outer-roll-gimbal error, also Ixsted in Table Bl. 

Star Tracker Errors 


The star tracker error consxsts of two basxc components. The noise equivalent 
error is random, being about 20 s^ (lo ) . The second error xs a function of the 
position of the star xmage xn the star tracker field of view. The star tracker 
procurement specification stxpulates that the star tracker shall measure the sensed 
star LOS with a total error not exceeding 60 (1 a) over the entxre fxeld of 

vxew. The per-axls requirement, normal to the tracker's boresight axis would 
therefore be 60/ , or about 42 s^. 

Although the star tracker is believed to perform somewhat better than the 
specxficatlon value, the 42 sec value will be used in this analysis. Discussion 
with W. Swingle, NASA/JSC, indxcated that star tracker accuracy Improves when the 
star sightings take place near the center of the star tracker field of view. The 
star tracker field xs a 10 x 10 degree square, the center of whxch being the 
star tracker's boresight axis. If star sxghtings are restricted to the central 
4x4 degree field, star tracker measurement errors are halved. 

IMU Gimbal and Resolver Errors 


These error source values are based on reference 5 data. In general the 
effect of these errors is a fxinction of gimbal angle. 

Referring to Page 7, where the pitch resolver angle is zero, the pitch to 
outer roll gimbal non-orthogonality error of 30 sec can be visualized as a slxght 
rotation of the pxtch gimbal about the stable member azimuth axxs. This introduces 
an error in the IMU azimuth readings. For non-zero pitch gxmbal angles, this non- 
orthogonality also affects the outer roll gimbal resolver readings. 



The resolver readings themselves are subject to several errors. The offset 
error is the error existing when the gimfial resolver reading should he zero. The 
resolver readout errors, over the range of glmbal angles readings (0 to 360 degrees) 
can be expressed as a random bias and a combination of sinusoidal harmonics whose 
magnitude and phase are random variables. The la magnitudes are given in Table B2. 
The phases can be considered uniformly distributed between 0 and 360 degrees. 

Resolver readout quantization is 20 sec. The maximum quantization error is 
therefore 10 sec. Treating the quantization error as imlformly distributed leads 
to a lo quantization error of 10 s^/ or about 6 s^. 

The ginbal and resolver errors affect the DHJ alignment systematically as 
a function of the ginibal angles. If all gimbal angles changed significantly 
between two readings, the readout errors might be considered independent random 
errors. If the change were small, then the readout errors would tend to behave 
as biases. In general, when the Orbiter changes its attitude, some of the gimbals 
might change significantly while other ginbal angles might change very little. 

Thus, the nature of these errors varies somewhat unpredictably. However, since 
both the in-orbit alignment and alignment transfer maneuvers involve large changes 
in Orbiter attitude, for this analysis the gimbal and resolver errors will be 
assumed independent random errors. 

Effect of Timing Errors 

Table B3 lists timing errors associated with the in-orbit alignment of the 
IMU. Only three error sources exceed 1 millisecond (ms). Of these, the star 
tracker data staleness error ranging from 10 ms to 52 ms predominates. The in- 
orbit alignment equations mechanized in the flight computer average 21 data points 
of star tracker and IMU gimbal angle data taken 160 ms apart. Thus, the averaged 
staleness error is about 30 ms. 


Timing errors affect the alignment only if the Orbiter is rotating durii^ the 
star sightings. If we assume, for example, the Orbiter is rotating 0.1 deg/sec 
during the measurements, an average error of .003 degrees or 11 s^ could be 
induced in the alignment. Presximably, the Orbiter 's angular velocity during the 
star measurements will be low, on the order of 0.1 deg/s^ or less. Hence, this 
analysis will ignore the effects of equipment timing errors. 
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TABLE B 3 


OB-ORBIT IMU ALIGNMENT TIMING ERRORS 


ERROR SOURCE 


MAGNITUDE COMFENSATABLE 


GPC time tag error 

IMU data input period skew 

IMU/star tracker data input 
deviation 

Star tracker data input period 
skew 

IMU 8x resolver data staleness 

IMU 8X resolver data staleness 
uncertainty 

IMU 8x resolver data skew 

IMU IX resolver data skew 

Star tracker data staleness 
uncertainty 


< U ms 

< 16.5 ps 

1782 ys 

< 16.5 ps 

5.608 ms 

< 100 ps 

< 100 ps 

< 12 ms 

10 to -52 ms 


Partially 

Bo 

Yes 

No 

Yes 

No 

No 

No^ (No 

No 


error effect) 




Since only the three most significant digits of these resolvers 
determine the octant of the IMU gimhal angle, no compensation is 


are used to 
necessary. 


OP POOR QUAUT^i " 
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APPENDIX C 


ORBITER IMU ALIGNMENT PROCESS 
AND 


ERROR ANALYSIS 






ALIGNMENT PROCESS 


The purpose of the alxgnment is to determine the transformation Tp^ from 
the stable member (present cluster) frame to the M50 inertial reference frame, 
obtain T^O, sightings are made on two different stars with a star tracker. The 
star tracker measurements, horizontal and vertical deflection angles relative to 
the star tracker boresight axis (body-fixed) , are coiribined with the IMU gimbal 
angles and known geometrical orientations between the star tracker and the IMU 
to compute unit star vectors in the IMU stable member coordinate system. Unit 
star vectors are available in the M50 coordinate system via a star catalogue. 
Therefore , let 

_r, be star vectors in M50 (published data) 

m, n, be measured star vectors (same stars as ^) 
in IMU stable member coordinates. 



s ^ 5 X j: 

3iar^Z 


Figure 2 Star LOS’s 


OF Pt)OR 



To obtain T pQ^, compute as follows: 


= unit (r^ X 

i X ^ 

form orthonormal triad in M50 


= unit (m X n) 
n^’ = ^ X m 

m, n^' , ^ form orthonormal triad in IMU stable member 
coordinates . 


Let matrix S = , where the indicated vectors form the columns of the matrix. 

S IS an orthogonal matrix (rotation transformation) . 

M50 

nijn' ,^| . M is an orthogonal matrix, T obviously transforms m 
to r , n* to s_’ , and to Therefore: 




M50 


T IS found by solving the matrix equation above. 

JrCr 

where is the transpose of M. For orthogonal matrices, the transpose is also the 
inverse transformation. 


ORIG^AL PAGE IS 
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RELATION OE THE ALIGNMENT ERRORS TO THE MEASUREMENT ERRORS 


The errors affecting the alignment accuracy are assumed to be due to the 
measurements by the star tracker and IMD and the mechanical alignment (geometrxcal 
orientation) ^of the star tracker to the IMU, Star catalogue errors are assumed 
neglxgible. Consider the figure below: 



Fiaure 3 Measurement Errors 

where ; 


per-axxs measurement errors e - , s . , are in the inertxally fixed plane formed 
by m and e - ti 2 > perpend^ular to the plane. A xs the angle between m 

and ^ and ^ represent the vector pointing errors xn m and _n respectxvely . 

The star LOS vectors contaxnxng the measurement errors are m + ^ and n + Sn . 
These unit vectors generate an orthonormal trxad of vectors sxmilar to m, 





The two triads, almost coin- 


This second triad we denote m + n'+ SnJ + SS. . 

cident, are related to first order as foTTows: 


r • ' 

1 1 

1 1 


1 ■%' 


r ' ' 1 

' ' I 

, I 


— 



bi n'|i 

' 1 

1 1 

J ) 


^ \ 

0.^ 

t 


' I 

I 

1 

1 


50^, 60^, 50£ represent small angular rotations about m, n_' , and respectively. 


The angular error vector 60 = Orbit er alignment error. 

We see that the matrix error equation is * 


1 1 

! 1 

1 1 


r 

o 

i 


~ ' 1 T 

: ; I 

t 1 

£m 1 1 

1 

. — 

£ 91 , 0 


B; n'li 

1 ] 

‘ 1 

1 1 

1 _ 



0 

\ 1 

1 1 

1 1 


The accuracy will be analyzed in the coordinate frame determined by the 
orthonormal triad m, n ' , U Thus , 


m = ( 1, 0, O) 
n = (Cos A, Sin A, O) 
n* = (0, 1, 0) 
i. = ( 0, 0, 1) 

^ = ( 0, -e^^) 


ORIGINAL PAGE 
OP Poor qxjautz 


— “ ~^nl ^nl "^n2^ 


ORIGINAL If ^ 
OF POOR 'QUAWy 
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¥e need to determine 6^ and SP . 

= nnit (m X n) = g " • £ (m x n) 

5i = - Cot A ( m X _n ) + Csc A ( ^ x n + m x ^ ) 

Substituting the vector components, performing the vector products, and combining 
terms give : 

= ( e s ^ Csc A - e _ Cot A, 0 ) 

■ — m2 n2 m2 

Next, we have 

n ' = X m 

6n * = 6J_ X m + X &m 

Substituting the vector components, performing the vector products and combining 
terms give: 


6n' = ( -e ^ , 0, -So Csc A + e ^ Cot A) 

— ml ’ ’ n2 'm2 

Substituting 6m , 5n * , 6i and m, n^’ , and in the matrix error equation above yields 


0 

~ S 



0 



0 

‘ 


— 


0 



0 





I 
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Thus : 



I 




s^yCotA - e^^CscA 

80 


= 

^m2 




^ml 


It is seen that 60jj^»and 60^ depend respectively on eJJ ^2 and emi These per- 
axis alignment errors exactly equal the per-axls measurement errors indicated. 

Small values of A degrade only the alignment error 60jjj about the m axis. The 
alignment errors 60jjt and 60^ about the other two reference axes n' and ^ do 
not depend on A. As A approaches 90 degrees 60ni approaches -en2> so 60m tends 
also to exactly equal a per-axis measurement error when the subtended angle is 
large. Thus, with A near 90 degrees, the Orbiter alignment errors are essentially 
equivalent to the measurement errors, which, intuitively, should not come as a 
surprise. 

Note the effect on alignment accuracy. This is because the "first" 

star m was chosen as one of the orthonormal triad vectors, while the "second" star 
n was not, being replaced by ii’ . 

In short, it is seen that the per-axis alignment error, for angle A in the 
vicinity of 90 degrees, essentially equals the per-axis measurement error. For 
Orbiter in-orbit alignment, the rule of thumb whereby the alignment error is 
computed as "^^times the measurement error is not valid. 

Alignment Error Statistics 

It will be seen later in this report, that the inertial per-axis measurement 
errors e-ml , e-m?. and eji 2 3^^®sums of errors stemming from the star trackers, the 
IMU’s and the geometrical misalignments of the IMU’s relative to the star trackers. 
We will assume Ejjj 2 , and normally distributed (central limit 

theorem) and, based on the data at hand, independent of each other. * 

It will also be seen that the geometrical misalignments essentially are* body 
fixed biases, during the period of an in-orbit alignment. Since e mi , . and £^2 

are referenced to an inertial frame, the question of Orbiter attitude maneuver 
policy arises. If the Orbiter performs standardized orientations in inertial 
space when sighting on given stars, then the body fixed biases will map the same way 
relative to inertial space. Hence, eml, em2> £n2 would have significant bias 
components for a given alignment procedure. We will, therefore, assume the 
measurement errors normally distributed about a bias, in general. 
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Denoting the error random component as r with standard deviation a and the 
hias as b, we have: 


ml 

^ml ^ml 

■m2 

''in2‘^^m2 

n2 = 

^n2 + ^n2 


Hence, the mean square errors are: 



•u 

~ r 

luj. 

“f" 

<X^ 

btyti 

<2^ 

•s=^ cr\ ^ 
ml 

"f- 

b 

mi 

X 

X 



4. 

4- 


m<i. 

— 

4- 


- 

~ 





! ^ 

X 


IX 

Q 

~ C 

+- 

[) 

r; (JN 

t- 


nn~ 

ni 

^U2- 

nx 


Also* 


^ ml hi2- 


= h b 

■* hfl fM2- 


— bnit b 


n<i- 


¥e now compute the mean square per-axis alignment errors. The result is. 


% — + £„IC£c"A 
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Substituting the measurement error 


statistics yields: 


^ ^ 



z 

^ -t 
mt 



Computing the error covariance matrix 



yields; 


Cov 



L 





0 


1 


'mx 


2, 

O 

hii 


0 


0 


O 


2 - 

(5N 

mi 


The variance of the total alignment error is the trace of the above matrix. 


The data on hand does not suggest significantly different statistics for 

E e Thus, we will assume the following: 
m2 ii2 



•2_ 

7_ 

A 

— 

^ — 
m2. 

^ ha- 

— 


ON ir 

ON 

A 

hil 


no. 


2. 

, 2, 

1 ^ 


^ >1 ” 

L 

TMO- 
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Hence ; 


0^ = 0 ^(Cot^A + Csc^A) + b^(+ Cot A + CscA)^ 

m — — 

0j = a ^ 


2 2 ,2 
0 = a + h 


In general, A ^ 90 degrees, so the alignment error about m is different from 
the alxgnment errors about n^' and For convenience, we now define the "averaged" 
per-axis alignment error to be 1/Y~3” times the total alignment error. Thus, the 
root mean square (RMS) per-axis alignment error is defined as 


per-axis SMS error = 




0^(1 + 2Csc^A) 


+ b 


[ 


2 + (+CotA+CscA) 



The standard deviation of the per-axis alignment error is; 

j 4 - 1 + 2Csc^A \ 

per-axis std. deviation “ 0 ^ ^ J 




As A approaches 90 degrees, the SMS per-axis alignment error approaches 


2 2 

0 + b , the SMS per-axis measurement error. 
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APPENDIX D 


OSBITER SENSOR BODY FIXED 


BIAS ERROR REMOVAL 






BODY-FIXED BIAS KEMOVAL 


Because each star tracker boresight axis is fixed relative to the Orbxter and 
because the star tracker fxeld of view is relatively narrow (horxzontal and vertical 
maximum deflections are i 5 degrees) , then a simple procedure can be invoked to 
remove the body-fixed alignment errors, as follows. 

After taking a star sighting simply rotate the Orbxter 180 degrees around the 
star LOS. Take a second readxng (same star tracker 1). For each reading, conqiute 
the star LOS unxt vector in stable member coordxnates. Average the two vectors by 
adding them together and divxding by 2. The effect of the body-fixed bias errors 
to first order is removed by the averaging process. 


Thxs may be clearly seen by considering the fxgure below, which represents 
the region of the celestial sphere in the star tracker field of view. Assume 
star 1 is in vxew. The unit vectors m, m^, and ^2 are directed from the viewer 
toward the field of view (into the page) . 



Figure Dl Bias Removal 


fxrst measurement, offset from the true value m by the 
body-fiiid"bxS\ and random error |2 |^“he biaf elfectroJ S2^are 

“f “Srrn: if r S - ..a. e„e. - , 

and random error ^2* averaged measurement is 


m 

— av 
m 

— av 
m 

— av 


%(m + b + 6r^ + m - b + ^ 2 ^ 
m + ^(drj^ + 
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,OP POOR QUAU3^ 



It xs seen that the bxas effect ^ is removed from The per-axis standard 

deviatxon of varies from 0 to a (where a is the per-axxs deviation of both 
6rj ^ and dr? ) , depending on the correlation between dr^ and dr? . 


For uncorrelated dr^ and dr? the per-axis standard deviation of m^^ is o/ViT 


The bias effect ^ could be estimated, if desired, by dxfferencing measurements 
m^^ and ^2 and dividing by 2. 


Si - 


—est 


H.2 - B + “ B k “ ^^2 

= _ m2> = ;b + ig (dr^^ - ^ 2 ) 


The per-axxs error xn would vary from 0 to a, dependxng on the correla- 

tion. For uncorrelated dr^ and dr? , the per-axis deviation xn the estxmate 
would be a/ \fz . 

Note that the 180 degree rotation need not be a precision maneuver. The rota- 
tion need not be strictly around the star LOS . It could be nominally around the 
tracker boresxght axxs, xf desired. It xs not necessary that the measurement after 
the rotation be taken with the star image in the same location in the tracker field 
of vxew as the fxrst measurement. It xs not important that the star be maintained 
xn the tracker fxeld of view durxng the maneuver. Also, rotation wxthxn il 3 degrees 
of the nomxnal 180 degrees is entxrely adequate, for practxcal purposes. Sxnce 
the star trackers poxnt essentxally along the Orbiter’s pxtch and yaw axes 
respectively, the Orbiter would be rotatxng about princxpal axes of inertxa during 
the rotations, rendering these 180 degree maneuvers relatively easy to perform. 

It should also be noted that the procedure above does not necessarxly remove 
angular bias errors about an axis coinciding with the star tracker boresxght axis. 
However, angular error effects about the boresight axis are reduced by at least 
an order of magnitude since the star tracker measurements (horizontal, vertical 
deflectxons) are restricted to + 5 degrees of the boresight axxs. Hence, to first 
order, the angular bxas error about the boresight axis has little significant 
xmpact on the Orbiter alignment accuracy. 
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ORBITER ALIGNMENT ACCURACY 


Two sets of results are presented. 

The first set of results, presented in Figure 5, pertains to the general 
case , where all error sources are treated as random with zero means . 

The second set of results, presented in Figure 6, pertains to a single in-orbit 
alignment process with the geometrical misalignment biases removed by averaging two 
measurements per star, measurements taken before and after rotating the Orblter 180 
degrees aromd each star LOS. The remaining random errors are considered uncorrela- 
ted; hence the averaged measurement deviation is a/^^~2 (68 sec/'^^Y'= 48 sec). 

For both sets, the effect of Orblter IMO gyro drift (.035 deg/hr, la) is added 
(in quadrature) to show how the in-orbit alignment degrades with time after it is 
initially achieved. Figures 5 and 6 are contours of error values, as functions of 
subtended star angle A and elapsed time after alignment takes place. Initial 
alignment error, as a function of star angle A, is indicated along the vertical 
axes, where the contours intercept the vertical axes at elapsed time equal zero. 

The allowable elapsed time, during which the Orblter ’s per-axis attitude error is 
less than 6 min (3 a) , is determined by moving horizontally to the right (from the 
initial alignment error at time zero) until the 6 min contour is intercepted. 

Figure 5 shows that with the misalignment biases present in the star measure- 
ments, a maximum time of about 18 minutes is available (A = 90 deg) before the 
Orblter attitude error exceeds 6 min. The star angle A must be x<rithin 90 "t 20 
degrees in order that the initial alignment error not exceed 6 min. 

Figure 6 shows that with, the jqis-slignment biases removed, a maximum time of 
over 52 minutes- is- available (A ?! 90 degl before the OrBxter attitude error exceeds 
6 min. The range of allowable s-tar angles is^ aften to be quite broad C9b 70 Deg.). 
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